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SYNOPSIS 



We study decoherence of one, two, and n non-interacting qubits. De- 
coherence, measured in terms of purity, is calculated in linear response 
approximation, making use of the spectator configuration. Monte Carlo 
simulations illustrate the validity of this approximation and of its exten- 
sion by exponentiation. Initially, the environment and its interaction with 
the qubits are modelled by random matrices. Purity decay of entangled 
and product states are qualitatively similar though for the latter case it is 
slower. 

For two qubits, numerical studies reveal a one to one correspondence 
between its decoherence and its internal entanglement decay. For strong 
and intermediate coupling to the environment this correspondence agrees 
with the one for Werner states, for initial Bell pairs. Using this relation we 
are able to give a formula for concurrence decay. In the limit of a large 
environment the evolution induces a unital channel in the two qubits, 
providing a partial explanation for the relation above. 

Using a kicked Ising spin network, we study the exact evolution of two 
non-interacting qubits in the presence of a spin bath. Dynamics of this 
model range from integrable to chaotic and we can handle numerics for 
a large number of qubits. We find that the entanglement (as measured by 
concurrence) of the two qubits has a close relation to the purity of the pair, 
and closely follows an analytic relation derived for Werner states. As a 
collateral result we find that an integrable environment causes quadratic 
decay of concurrence as well as of purity, while a chaotic environment 
causes linear decay. Both quantities display recurrences in an integrable 
environment. Good agreement with the results found using random ma- 
trix theory is obtained. 

Finally, we analyze decoherence of a quantum register in the absence of 
non-local operations i.e. n non-interacting qubits coupled to an environ- 
ment. The problem is solved in terms of a sum rule which implies linear 
scaling in the number of qubits. Each term involves a single qubit and its 
entanglement with the remaining ones. Two conditions are essential: first 
decoherence must be small and second the coupling of different qubits 
must be uncorrelated in the interaction picture. We apply the result to 
the random matrix model, and illustrate its reach considering a GHZ 
state coupled to a spin bath. 

PACS numbers: 03.65.Yz, 03.65.-w, 03.65.Ud, 05.40.-a 

Keywords: entanglement, random matrix theory, purity, decoherence, 

concurrence, quantum memory, quantum register, GHZ. 



RESUMEN 



La decoherencia de 1, 2 y n qubits no interactuantes y posiblemente en- 
lazados, medida en terminos de la pureza, es calculada usando respuesta 
lineal y el concepto de configuracion de espectador. A traves de simu- 
laciones de Monte Carlo exploramos la validez de la aproximacion y su 
extension mediante exponenciacion. Inicialmente, modelamos la inter- 
accion y el medio ambiente con matrices aleatorias (MA). 

Para 2 qubits, en el modelo MA, el enlazamiento interno y la decoheren- 
cia tienen una relacion uno a uno. Esta relacion, para acoplamientos 
moderados y fuertes, coincide con la relacion correspondiente para esta- 
dos de Werner, si la condicion inicial es im par de Bell. Mediante esta, 
obtenemos una formula explicita para el decaimiento del enlazamiento 
interno. 

Introducimos un modelo de Ising pateado (MIP) para estudiar un grupo 
de espines, acoplado debilmente a un bano de espines. Este modelo pre- 
senta dinamicas integrable, mixta y caotica para diferentes parametros. 
Observamos nuevamente la relacion entre la decoherencia y el enlaza- 
miento interno, obtenida con el modelo MA. Inicialmente, tanto el en- 
lazamiento como la decoherencia decrecen cuadraticamente/linealmente 
para el caso integrable /caotico. En el caso integrable, si el acoplamiento 
es a im extremo de una cadena, ambas cantidades presentan compor- 
tamientos periodicos. Comparamos cuantitativamente el decaimiento de 
pureza en el modelo MA con nuestro sistema dinamico. Los resultados 
positivos demuestran la validez del modelo estocastico. 

Finalmente, analizamos una memoria cuantica expuesta a decoheren- 
cia. Nuestro resultado, que asume altas purezas e independencia de 
los acoplamientos en la imagen de interaccion, permite expresar la de- 
coherencia como una suma de terminos que involucran un solo qubit 
y su enlazamiento con el resto de la memoria. Aplicamos el resultado 
al modelo MA, generalizando nuestros hallazgos. En el MIP, observa- 
mos que incluso en situaciones integrables, los requerimientos se pueden 
cumplir. 

Numeros PACS: 03.65.Yz, 03.65.-w, 03.65.Ud, 05.40.-a 
Palabras clave: enlazamiento, enmaranamiento, RMT, pureza, decoheren- 
cia, concurrencia, memoria cuantica, GHZ, Bell. 
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Chapter 1 

Introduction and fundamental 
tools 

Studying decoherence of one, two, and n qubits has a wide scope of appli- 
cations due to the huge interest in implementing "quantum technology". 
The limiting factor for building this technology is the sensitivity of quan- 
tum systems to undesired perturbations/ coupling. Moreover, coupling 
to external degrees of freedom is fundamentally inevitable. Understand- 
ing its behavior is crucial to tame its effects. 

Since a long time the coupling of quantum systems to external degrees of 
freedom has been studied (see e.g. llvN55l IEve57l IAlb92l lArb 93l). Philo- 
sophical aspects of quantum mechanics (the measurement process and 
the emergence of the classical world) are deeply connected with the prob- 
lem. Some particular models designed for specific applications have been 
developed, but the favorite for general purposes (by far) is the Caldera- 
Leggett model PCL83I in which the external degrees of freedom are mod- 
eled by a set of harmonic oscillators. Good agreement with the experi- 
ment has been observed, e.g. IIKD98ilWFL02llZCP+07J . 

Regarding applications to quantum information some progress has been 
made. A big amount of literature exists and some aspects have been 
demonstrated experimentally. Most theoretical studies use Caldeira-Le- 
ggett like models; others explore the consequences of using/dropping 
the usual Markovian approximations. Some use very particular models 
either to obtain explicit analytical results or to study specific experimen- 
tal situations. A more general picture is thus desirable to gain a deeper 
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understanding of the physics governing decoherence of quantum infor- 
mation systems. A few comments on some relevant papers on the field 
are useful to have an idea of the situation in the literature. This list of 
papers is not meant to be complete, it is just intended to give a brief 
overview of what people are currently working on, in relation to this 
topic. 

• In a series of papers IIYE021 IYE031 IYE04| Yu and Eberly explore how 
two qubits (in particular their entanglement properties) are affected by 
broadband bosonic reservoirs that induce vacuum noise, phase noise, 
etc. They find that entanglement after a finite time goes identically to 
zero. 

• In IIBra 061 a more complicated study is done using again a traditional 
harmonic oscillator bath. There, the off diagonal elements of the re- 
duced density matrix (sometimes called "coherences") are analyzed. 
Interestingly, the author discovers that decoherence is determined by 
a generalized Hamming distance. He introduces some primitive spec- 



tator (see sec. 3.1 1. 



In ||Ged06ll the author uses a spin bath as an environment. He studies 
concurrence of some Bell pairs. The results, though interesting, are 
quite model dependent. In lLDK+051 the relation between integrabil- 



ity and decoherence is studied for a spin bath environment, very much 
in the spirit of our results IIPS06I . Studies of specific spin-bath envi- 
ronments, aiming to understand decoherence in experimental qubit 
realizations are lldSDS03blldSDS03all and IISLH+041 . 



In another interesting article IIGMCMB07I , Garcia-Mata et al. use 
semi-classical considerations to study multi-particle qubit entangle- 
ment. They analyze how entanglement is affected both by the phase 
space structure and by the kind of noise applied. 

In IIFFP04II the authors propose two measures for decoherence. These 
measures are additive, a property that is important to express the total 
decoherence in terms of the decoherence of each qubit. 

Some efforts to understand the implications of the Markovian approx- 
imation are made in I ALKH02,.LKA+04J , where it is shown that under 
certain circumstances that approximation leads to very inaccurate re- 
sults. 

In |MCKB05|I the authors use random Hamiltonians (though not with 
the minimum information properties of the classical ensembles l|Bal68l ) 
to study entanglement decay of n qubits. Using Markovian approxi- 
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mations, they arrive to time independent Linblad equations. They ob- 
tain multi-exponential decay. They are able to analyze the differences 
between W and GHZ states. 

• An isolated and possibly premature (due to the interests of the com- 
munity) study is worth mentioning. In l|MPK88| Mello et al. analyze 
the relaxation rates of a single 1/2 spin particle using a random matrix 
model. 

We are pioneering the use of Random Matrix Theory in the field of quan- 
tum information theory. Some previous work has influenced consider- 
ably our research namely [GPSQH IGS03I , where random matrices were 
used to analyze decoherence and fidelity in general quantum systems. 

Joseph Emerson has developed methods to create random unitary opera- 
tors (in the spirit of the CUE) using random gates (ironically a non-trivial 
task to implement efficiently) lEWS^OSj . He has also explored the utilities 
of such random unitary operators in fidelity and local density of states 
estimation IIELPC04I . Other uses of randomization in quantum protocols, 
like diminishing the effects of static perturbations, have been introduced 
in ||KAS05[[PZ01|| and further explored in [KA06|. Some people have also 
used random matrices to study fidelity decay I^GPSZ06l IFFS04| . 

This thesis is based mainly on four publications IIPS071 IPGS071 IPS06[ 
|GPS07|. We shall not follow the chronological order as the logical order 
will result in almost opposite. The reasons are clear. As you gain insight 
in the field, things become clearer, concepts become more elaborate and 
distilled and thus more suitable for understanding. 



On the structure of the thesis - 

In the remaining part of the introduction we explain some concepts and 
tools used along the thesis. Some material is not new, and is not appro- 
priate for publication in a journal as it only contains a review of known 
things. However, for the reader a coherent presentation is always handy. 
A main concept used and studied during this thesis is entanglement. It 
is understood here both as a resource (to perform quantum information 
tasks) and as a cause of decoherence. We shall define it and discuss the 
two ways that we understand it. Next we explain the spectator configura- 
tion, which is an original tool exploited during most of this work. Finally 
we introduce Random Matrix Theory (RMT). 
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We next proceed to analyze decoherence of quantum systems. We first 
explore the single qubit case (sec. |2j|. The detailed mathematical deriva- 
tion of the formulae used is given in appendix |A] Both the time reversal 
invariant (TRI) case and the non-TRI case are analyzed in detail. We then 
consider the two qubit case (sec. |3]). Different configurations correspond- 
ing to different physical situations are analyzed. Again TRI and non-TRI 
cases are also studied. The effect of internal entanglement proves im- 
portant and provides interesting effects. In sec. |4] the relation between 
decoherence and entanglement is studied. Sections |2| |3| and |4] are based 
on IIPGS07I , though the basic idea was introduced in IP SOTl . 

In sec. |5] we give an example of how some of the concepts can be applied 
to a simple model: the kicked Ising spin chain. Though this chapter is 
based on |PS06L major modifications have been introduced with the aim 
of getting closer to the RMT models. 

Finally in sec. |6] we use some of the results obtained during the thesis to 
analyze decoherence of an n qubit register. We shall use both the RMT 
model and the KI spin chain to discover and understand the reach of the 
results. This chapter is based on [GPS07|. 

In appendix |A] we perform the main RMT calculation, for a single qubit 
in the spectator configuration. In appendix |B] we explain how to imple- 
ment numerically the kicked Ising model in an efficient manner. The next 
appendix (jCj) explains a way of extending some analytical results via ex- 
ponentiation. This heuristic result is tested throughout the thesis for most 
Monte Carlo simulations. We then discuss some technical aspects of both 
entanglement (appendix |D| and random matrix theory (appendix |E]|. Re- 
garding entanglement, we discuss the definition of entanglement in more 
general systems than the ones discussed in this introduction, and the 
physical meaning of concurrence. For random matrix theory we mention 
the physical justification of the ensembles, relations among their matrix 
elements, and other formulae used during the thesis. In the last appendix 
we give the double integral of the form factor for a particular ensemble 
and a simple proof of the Born expansion for the echo operator. 

I do hope you enjoy and have a nice time with this piece of work. 
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1.1 Entanglement 

Though in mathematical terms entanglement is trivial to define (once 
the basic tools of quantum mechanics are introduced), its consequences 
challenge many deeply rooted (mis)conceptions about reality. Here we 
do not wish to discuss how the existence of entanglement affects our 
understanding of reality; this is a difficult topic outside the scope of this 
work, and even Einstein was puzzled by its consequences. We limit our 
selves to define and discuss briefly entanglement and how to quantify it. 

A pure state of a quantum composite system is said to be entangled when 
it is not the "sum" of its parts (technically we mean tensor product). 
To be more precise, let our Hilbert space Ti. be composed of two parts: 
H = Ha ® Hb- If, given a state \xp) & H, there exist \tpA) € Ha and 
|i/?b) e Hb such that 

\xp) = IxPa) ^ \^Pb) (1.1) 
it is said that \tp) is separable or unentangled. Conversely, if 

\XP) ^ \lpA)^\tpB),y{\lpA) ^Ha, \iPb)^Hb) (1.2) 

it is said that | xp) is entangled. In other words | xp) is entangled if and only 
if 

I (t/;| ^tr^|i/^)(t/;| strait/;) (t/;|. (1.3) 

It is quite easy to show the existence of entangled states. The simplest 
case can be constructed when dimHA = dimT^g = 2, i.e. when Ha and 
Hb represent qubits. Let {|0), |1)} be an orthonormal basis in each space. 
The Bell state 

IBell) = \'»^\'»^^}^m , „ ,1.,) 

is entangled. Assuming the existence of cx.a,^a,'X-b,^b € C, such that 
|Bell) = (^aIO) + ® (aslO) + |6b|1)) results in a contradiction. For 

multipartite mixed systems a generalization of the definition of entangle- 
ment is straightforward. See appendix [D| for details. 

In order to get deeper insight in the entanglement properties of pure bi- 
partite states it is convenient to use the Schmidt decomposition IISch07[ 
INCOO I. Given a state |t/j) in a bipartite space Ha ^ Hb, there exist or- 
thonormal states {|/a)} in Ha and {|zb)} in Hb such that 

min{ dim Ha ,dim Hb} 

W)= £ A;|za) ® |/b) (1.5) 

! = 1 
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and < A,- < 1, with I^, A? = 1. The numbers A, are called Schmidt 
coefficients and play an important roll in entanglement theory. Consider 
an orthonormal (and complete) basis that diagonalizes pA = trg \xp){\p\. 
We choose that basis to be jz'/i); its existence is guarantied by the spec- 
tral theorem. We can then write \tp) = Ylil^A) ^ Ih), but since pA = 
E; |^a)(mI(zb|/b) must be in fact diagonal, then (zb|/b) o« ^ij- Using some 
CK |zg) such that (zbI^b) = 1 arid suitably choosing its phases we can 
write eq. \1.5\ . The sharp reader will notice that the Schmidt coefficients 
are the square roots of the eigenvalues of the reduced density matrix of 
any of the two subsystems. 

The Schmidt coefficients are unique for each pure state. From the argu- 
mentation we can see that Pa and ps have the same eigenvalues (and with 
the same degeneracy) except for | dimH^ — dimT^el zeros. Determining 
whether a pure state is entangled or not is an easy task. From the previ- 
ous paragraph one can see that a state is not entangled if and only if one 
of the Schmidt coefficients is one (implying that the others are zero). 

1.1.1 Decoherence as entanglement 

Decoherence can be seen as entanglement with the environment |Zur03[ 
IZur91l . 

Quantum correlations [in our language, entanglement] can 
also disperse information throughout the degrees of freedom 
that are, in effect, inaccessible to the observer iZur91| . 

Though a big debate has been issued since the formulation of that para- 
digm, it is now generally accepted. 

We now discuss an example to explain the previous statement. To present 
the key idea it is enough to consider a central system, composed of a sin- 
gle qubit, and an environment alone; in the original formulation a mea- 
surement apparatus was also involved to allow the analysis of the "col- 
lapse" of the wave function after a measurement process. In this example 
the environment has three characteristics: large dimension, uncontrol- 
lable dynamics, and no possibility of being observed. We assume some 
interaction between the central system and the environment. Consider an 
initial state which is (i) separable with respect to the environment and (ii) 
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(1.6) 



where [0) and [1) form an orthonormal basis for the qubit, a and /3 are 
complex numbers, and \(p) is the initial state of the environment. Assume 
that the interaction depends on the state of the qubit, e. g. a controlled- Lf. 
After some time, due to the interaction, the state will be 



As the dimension of the environment is big, the states of the environment, 
after some time scale, will be approximately orthogonal: {(pi\(pj) ~ Sij. Of 
course this last statement is not fulfilled for an arbitrary interaction, but 
precisely the basis (regarding the qubit) in which this condition is fulfilled 
will determine the preferred basis which determines the pointer states. 

To quantify the degree of entanglement of a bipartite system, in a pure 
state, we make use of the Schmidt coefficients. Adding any convex func- 
tion of these coefficients is enough. We use the sum of their squares, as 
it induces a very simple formula (other common choice, instead of x^, is 
xlogx which induces the von Neumann entropy). This measure we call 
purity. Thus, for a given density matrix p, its purity is defined by 



This quantity is 1 for pure states {p = \'^p){^\), less than one for mixed 
states, and reaches a minimum of 1/N (where N is the dimension of p) 
for the completely mixed state 1/N. If the partial trace with respect to 
an environment is represented by tre, a measure of decoherence is then 
P{p = tre \ ^){''P\)- An important practical advantage of this measure is 
that one does not need to evaluate the Schmidt coefficients of the density 
matrix p. 

Other views of decoherence are common in the literature. Consider a 
qubit in an initial state \ip) = (|0) + \l))/^/2. Its corresponding density 
matrix is 



Assume we have pure dephasing (no amplitude damping). Tj^ically, 
what will happen is that the off diagonal elements will decay exponen- 
tially. That is, its time evolution will be 



|t/;(bigO) = «|0)|<^o)+|6|l)|<^'i). 



(1.7) 



P(p)=trp2. 



(1.8) 




(1.9) 




(1.10) 
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Inspired in this behavior one can relate decoherence to the norm of the 
off diagonal term. We define 



D(p) =41^21'. (1.11) 



For states of the form eq. (1.10 1 we obtain the formula D{p) = P{p). 
However, in the general case, information about one only gives partial 
information about the other; for an arbitrary one qubit density matrix, 
< D{p) < P{p), and thus one can have a completely pure state with 
D = 0. This quantity is used frequently as is easy to calculate and is 
related to the interference fringes shown in the very popular cat states in 
phase space, see e.g. IIZur03l page 742. A big disadvantage of using D is 
that it is a basis dependent quantity. Internal dynamics may produce a 
decay of the off diagonal elements of p and thus of D. Moreover, if one 
studies an (n > 2)-level system the situation becomes more complicated. 
Purity on the other hand works for a much wider class of systems. 



1.1.2 Entanglement as a resource 

It is not difficult to understand that entanglement is a (quantum) resource, 
since already classical correlations are an important (classical) resource 
used extensively in classical cryptography. Entanglement, as the quan- 
tum correlation, brings up richer possibilities. In general, controlled en- 
tanglement can be used for the following: 

• Teleportation: It is the most celebrated application, due to its spec- 



tacularity and simplicity [ BBC"^93| . The transfer of an unknown 



quantum state can be achieved using an entangled state, local oper- 
ations, and classical communication. 

Quantum computation: It is a controversial subject whether en- 
tanglement is essential for quantum computation, but so far it has 
been demonstrated that for an exponential speedup in pure state 
schemes, entanglement is necessary (see |JL03| ). 

Communication: Both quantum and classical communication can 
benefit from entanglement. In particular, quantum key distribution 
extensively uses this resource [NCOOJ. 

Quantum-Enhanced Metrology: It is shown that the signal /noise 
ratio can be increased qualitatively nGLM06l IGLM04I if one uses 
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entangled states. Thus the use of highly entangled states shall be 
mandatory for precise measurements. Generalizations of the ideas 
developed in this area can be used to build quantum positioning 
systems (in analogy to GPS), enhanced radars, and for clock syn- 
chronization. 

Still the field is quite young and ideas for exploiting entanglement are 
emerging at this point. New technology is arising: some quantum ran- 
dom number generators are already available as USB gadgets. In the near 
future many expected, and unexpected, technologies are going to be pro- 
posed and, no doubt, realized. Thus it is a major concern to be able to 
understand, quantify, and control internal entanglement. 

One of the first tasks of quantum information theory was to quantify the 
degree of entanglement. It was soon realized that, in general, this was 
a complex task. We now know, for example, that for general systems, 
entanglement induces only a partial ordering. 

We now focus on the simplest possible scenario that allows entanglement, 
a two qubit system. Four conditions must be fulfilled by an entanglement 
measure: (i) It must have a value between zero and one. It is zero for 
separable states and one for Bell states, (ii) Any local unitary operation 
leaves entanglement unchanged. This condition can be seen as an invari- 
ance of the measure under a local change of basis, (iii) Local operations 
plus classical communication cannot increase entanglement. I.e. to create 
entanglement we need genuine non-local quantum operations (say inter- 
action, skew measurements, etc), (iv) The entanglement measure must 
be a convex function. This condition says that entanglement will not in- 
crease when mixing ensembles. These four conditions can be generalized 
to more complex systems (multipartite or higher dimensional systems). 

Several measures of two qubit entanglement fulfill these conditions, how- 
ever these measures do not provide exactly the same ordering of states 
IIVADMOll . In this work we shall use the concurrence. The first reason be- 
ing that it is straightforward to compute. Some measures of 2 qubit mixed 
state entanglement require explicit maximization over high dimensional 
continuous sets. Though in the definition of concurrence (via the entan- 
glement of formation) a maximization is required, the problem is solved 
in a general fashion, and a closed formula is given. See appendix |D] for 
details. The second reason being that it is used widely in the community 
of quantum information, both by theoreticians and experimentalists. 
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Concurrence C of a two qubit density matrix p is 

C(|0) =max{0,Ai- A2- A3- A4} (1.12) 

where A, are the eigenvalues of the matrix -y^ pi^Cy ® (ry)p*{ay ® cry) in 

non-increasing order. The superscript * denotes complex conjugation 
in the computational basis and dy is a Pauli matrix. Furthermore, con- 
currence fulfills all conditions of a legitimate entanglement measure dis- 
cussed at the beginning of this section. 



1.2 The spectator configuration 

One of the important contributions of this work is the concept of spectator 
configuration. During the development of the thesis the concept was 
discovered, and its potential is exploited here. On one hand, it allows 
to enclose all the calculations in a single one, thus simplifying greatly 
the technical details. On the other, it enables to extend easily our results 
from one and two qubits to n qubits. Its full potential has not yet been 
exploited, but we hope that the commimity will take advantage of this 
concept. 

The concept involves the following Hilbert spaces: 

• The spectator space Hs- The spectator is not coupled to the other 
spaces. However, it can be correlated initially via entanglement with 
the interacting space. 

• The interacting space Hi. This subspace is dynamically coupled to 
the environment and is initially entangled with the spectator sub- 
space. 

• The environment space He- This space typically has a large di- 
mension (though this is not essential at this point). It is initially 
decoupled (in the sense of entanglement) to the rest of the system. 

• The central system He- It is the tensor product of the spectator and 
the interacting space: He = 'Hs'Si Hi. 



The whole Hilbert space H is the tensor product of all spaces, namely 

H = He®Hc = He^Hs^Hi. (1.13) 



1.2. The spectator configuration 
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Additional to the Hilbert space structure, the spectator configuration, as 
explained above, has a characteristic Hamiltonian: 

H = He O li,s + Hs O le,i + Hi ls,e + Is O We,i (1.14) 
= He + Hs + Hi + We,i. 

The indices indicate the spaces in which the operators act. Where there is 
no danger of confusion, the identities are dropped as in the last equality 
of eq. ( 1.14| . The first three terms represent local Hamiltonians in each 



proper subspace, and the last term is an interaction between the inter- 
acting subspace and the environment. The different parts of the Hamil- 
tonian need not to be time independent; in this work the parts devoted 
to random matrix theory deal with time independent Hamiltonians. The 
parts studying the kicked Ising model and the n-qubit chapter use a time 
dependent Hamiltonian example of ( 1.14[ l. 



The initial condition is always separable with respect to the environment: 

\xp{t = 0)) = \xPe) ^ \tpc), \tpe) e He, \tpc) G (1-15) 

but not necessarily with respect to the spectator space. If there is separa- 
bility between the interacting and the spectator spaces, the problem triv- 
ially separates and we can consider then 2 completely decoupled prob- 
lems, one in Hs and another in He ^ Hi. The condition of the environ- 



ment being pure in eq. (1.151 is technical. Some calculations have been 
done with mixed states in the environment(s) yielding similar results. At 
this point we could formulate, with no problem whatsoever, the model 
with a mixed environment, but since for further considerations it is con- 
venient to have a pure state we keep it that way. 

The Hamiltonians that are going to be analyzed during the thesis are 
not always of the form eq. ( 1.14| , but do have a particular structure due 



to the structure of the underlying Hilbert space. This structure is of the 
form given in eq. ( |1.13 1, but with the interacting space being composed of 



different independent non-interacting groups. For this configuration, the 
environment will be coupled to all groups independently. Under some 
general conditions we shall be able to decouple this complicated problem 
into many spectator problems. 

Two explicit examples of a simplification of the problem using the spec- 
tator configuration are given when we have 2 qubits as the central system 
in the context of random matrix theory (sec. |3.3| and sec. |3.4| . Section |6] 
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separates the decoherence problem in spectator configurations in a gen- 
eral fashion, and exemplifies the results with both random matrix theory 
and the kicked Ising spin chain. As the reader can notice, we shall follow 
a line of argumentation that will build step by step the general case. Dur- 
ing the thesis we shall first study the simplest case (one qubit), then the 
next in complexity (two qubits) and finally explain a possible way to use 
the tools developed in a more general way. 



1.3 Random matrix theory: a tool 

All I know is I know nothing. 
Socrates 



Some people in the field consider the start of Random Matrix Theory 
as being the paper by John Wishart [,Wis28 l (who, incidentally, died in 
Acapulco). There he introduces random matrices with invariant measure 
under basis transformation. His objective was to analyze multivariate 
data. Others say that the landmark was placed by Elie Cartan in an old 
(and often forgotten) paper |Car35] . He introduces explicitly the circu- 
lar ensembles, with invariance properties with respect to (usually) group 
operations, to generalize the integral theorem of Cauchy. Mehta IIMeh91ll 
was the first to calculate many of the mathematical properties of the clas- 
sical ensembles [CarSSJ. 

There was not much development of RMT outside mathematics, until 
Wigner published his famous papers ||Wig51| |Wig55| pioneering its use 



in physics. These papers contain two important aspects. The first one is 
the idea to study statistical properties of the resonances of complex nuclei 
instead of studying its particular properties. This is in perfect analogy 
with statistical mechanics: One does not care about the particular posi- 
tion of the system in its phase space but rather about its thermodynamical 
(statistical) properties. The other idea was to use an ensemble of matri- 
ces to describe the system. This is again done in analogy with statistics 
mechanical, but it is conceptually quite different. Instead of performing 
averages over the phase space, one does averages over the space of sys- 
tems. Wigner was successful in describing some experimental findings, 
and later evidence showed the wide scope of applicability of his idea 
llBFF+81 1 IGMGW981 JPAl . 



1.3. Random matrix theory: a tool 
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Figure 1.1: Cartoon of an idea exploited in this thesis. The statistical 
properties of particular system (on the left, say a complicated 
billiard) are studied using the properties of an ensemble of ran- 
dom matrices (on the right). 



A next revolutionary step was marked by the papers by Casati et al. 
IICGVG80 I and Bohigas et al. IIBGS84 1. There they conjectured that quan- 
tum systems whose classical counterparts are chaotic have a spectrum 
whose fluctuations resemble those of the appropriate classical ensembles. 
The revolutionary aspect of this conjecture is that it does not require a 
complex composition of the system [i.e. many bodies), but only com- 
plexity in its dynamics. Vast numerical evidence favoring this conjecture 
is available IIGMGW98[ [JPA| , but a precise understanding {i.e. a globally 
accepted proof) is yet outstanding. What about systems with no clas- 
sical correspondence? Defining chaoticity in this case is cumbersome. 
We shall keep an oversimplified definition: quantum chaotic systems are 
those that exhibit fluctuations in its spectrum similar to those observed in 
the appropriate RMT ensemble. Alternatively one could say that quan- 
tum chaotic systems are those for which the correlations of most pair of 
observables decay to zero at large times. 



In recent experiments IISKK+OOi IZCZ+04[ lHSKH+05[ lHHR+05 1. it has 
been demonstrated that it is possible to protect ever larger entangled 
quantum systems, often arrays of qubits, ever more efficiently from deco- 
herence. A close connection between the dynamics of fidelity decay and 
decoherence has been shown in some instances [ CPW021 |KJZ02[ ICDPZ031 
IGPSS04L which suggests to apply methods successful in one field to t he 
other. In that context, a random matrix description ||GPS04l lGPSZ06t is 
accessible and very effective in describing experiments |SSGS05l ISGSSOSi 
IGSW06I . Based on this success of random matrix theory, we shall use it 
to model decoherence |GPSZ06[ [GSOZJ of qubit systems I PS071 , assuming 
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complicated dynamics in the environment, and a complicated coupling 
(in the interaction picture). 

Three perspectives make such a random matrix treatment particularly 
attractive. First, reduction of decoherence may, in some instances, be 
achieved by isolating some "far" environment (including spontaneous 
decay) to a degree that it can, to first approximation, be neglected. Then it 
can happen that the Heisenberg time of the relevant "near " environment 
is finite on the time scale of decoherence. In such a case it becomes 
relevant that RMT shows, in linear response approximation, a transition 
from linear to quadratic decay at times of the order of the Heisenberg 
time. This behavior is seen with spin chain environments [see sec. |5], and 
is essential for the success of the theory in describing the above mentioned 
experiments of fidelity decay. Note also that the concept of a two stage 
environment has been used for basic considerations IIZur03l . Second, the 
long term goal must be to describe in one theory the decay of fidelity 
that includes undesirable deviations of the internal Hamiltonian of the 
central system (already done), together with decoherence (done in this 
thesis). Third, random matrix descriptions include some aspects of chaos 
or mixing that are essential in the above experiments and may be useful 
for application to quantum computing [PZOl . F FS04 , CPZ05J. 

For some technical aspects of random matrix theory, including Heisen- 
berg time, GOE and GUE ensembles, form factor and density of states, 
please go to appendix |E] 



Chapter 2 

One qubit decoherence 



In this chapter we analyze decoherence of a single qubit. We focus on 
weak coupling of the qubit to an environment. We shall use the cor- 
relation function approach proposed for purity decay in echo-dynamics 
IIPS02I , treating the coupling as the perturbation. The linear response 
approximation will be sufficient. In this approximation the ensemble av- 
erages, which we have to take in any RMT model, are feasible though 
somewhat tedious. Exact solutions, which exist in some instances for the 
decay of the fidelity amplitude IISS051 IGKP^oH, seem to be out of reach 
at present, because they require the evaluation of four-point functions. 

The general program is as follows. Assume that the qubit is initially in 
a pure state, and evolves under its own local Hamiltonian. The qubit is 
coupled via a random matrix to a large environment in turn described 
by another random matrix. The coupling to the environment gives rise to 
decoherence. Averaging both the coupling and the environment Hamilto- 
nian over the RMT ensembles yields the generic behavior of decoherence 
of the qubit. 

We present a detailed analysis using both the Gaussian unitary (GUE) 
and the Gaussian orthogonal (GOE) ensembles IICar35[ |Meh91| for the 
description of the environment and the coupling. The two ensembles 
correspond to time reversal invariance (TRI) breaking and conserving dy- 
namics respectively. 



In sec. 2.1 we shall state the model, recall the linear response formalism 



for echo dynamics, and show how it can be adapted to forward evolu- 
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tion. In sec. 2.2 we discuss how to express the problem in terms of echo 
dynamics. In sec. 2.3 we give the general solution, arising from the calcu- 
lations done in the appendix |A] The analysis for the GUE case is given in 



2.4 whereas the one for the GOE is given in 2.5 



2.1 The model 



We describe decoherence by considering explicitly additional degrees of 
freedom (henceforth called "environment") which are interacting with 
the qubit. The Hilbert space studied in this section is 

n = n^(^ He, (2.1) 

where Ti.i (of dimension two) and Tig (of dimension Ng) denote the Hilbert 
spaces of the qubit and the environment, respectively. The Hamiltonian 
is of the following form 

Ha = Hi (g) le + li O He + Wi^e (2-2) 
= Hi + He + AVi^e ■ (2.3) 

Here, Hi represents the Hamiltonian acting on the qubit. He the Hamil- 
tonian of the environment, and Vi^e the coupling between the qubit and 
the environment. Notice how the indices in the operators indicate the 
spaces in which they act. The real parameter A controls the strength of 
the coupling. We shall study the time evolution of an initially pure and 
separable state 

\xp{t = 0)) = \iPi) \tpe) , (2.4) 

where \xpi) ETii and |t/^e) S 'He- At any time t, the state of the whole sys- 
tem is thus \ip{t)) = exp{—itHx)\Tp{0)), and the state of the single qubit 
is tre (t/;(f) |. As time evolves, the qubit and the environment get 

entangled, which means that after tracing out the environmental degrees 
of freedom, the state of the qubit becomes mixed. 

At this point we wish to compare this model with the spectator model. 
We can arrive to the one studied in this chapter from two different direc- 
tions. One is if we eliminate the spectator. The other is if we consider a 
separable (with respect to the spectator) situation [i.e. if in eq. ( 1.15| we 

let |i/;c) = \xps) ^ \tpi)l 

We describe both the coupling and the dynamics in the environment 
within random matrix theory. To this end, Hg and Ve,! are chosen both 



2.2. Echo dynamics and linear response theory 
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from either the GUE or the GOE, depending on whether we wish to de- 
scribe a TRI breaking or TRI conserving situation. The Hamiltonian H-[ 
implies another free parameter of the model, namely the level splitting 
A of the two level system representing the qubit. While the state of the 
qubit \tpi) implies more free parameters in our model, we assume the 
state of the environment {xpe) to be random. This means that the state is 
chosen from an ensemble which is invariant under unitary transforma- 
tions, and is fully consistent with our minimum information assumption. 
In practice, this means that the coefficients are chosen as complex random 
Gaussian variables, and subsequently the state is normalized. 



2.2 Echo dynamics and linear response theory 

We shall calculate the value of purity as a function of time analytically, 
in a perturbative approximation. As we want to use the tools developed 
in the appendix |A] for a linear response formalism in echo dynamics, we 
must state the problem in this language. To perform this task it is useful 
to consider the above Hamiltonian [eq. ( |2.2| ], as composed by an unper- 
turbed part Ho and a perturbation W. The unperturbed part corresponds 
to the operators that act on each individual subspace alone whereas the 
perturbation corresponds to the coupling among the different subspaces; 
i.e. Ho = He + Hi and V = Ve,i- 

We write the Hamiltonian as 

Ha = Ho + AV, (2.5) 
and introduce the evolution operator and the echo operator defined by 

U;,{t) = e-''''^', Ma(0 = ^o(0'^a(0. (2.6) 

respectively {h = 1 during all the thesis). The echo operator receives its 
name because it evolves a state forward in time with a perturbed operator 
and backwards with an unperturbed one. For the calculation of purity 
at a given time t, we replace the forward evolution operator Ha by the 
corresponding echo operator Ma- Even though the resulting states are 
different, i.e. 

pit) = tre,e' U;,{t)pUl{t) ^ p^{t) = tre,e' M;,{t)pMl{t), (2.7) 
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they are still related by the local (in the qubit and the environment) uni- 
tary transformation Uo{t). Since local transformations do not change the 
entanglement properties, it holds (exactly!) that 

P{t)=P[p{t)]=P[p''{t)]. (2.8) 

This step is crucial, since the echo operator admits a series expansion with 
much larger range of validity (both, in time and perturbation strength). 
However the numerical simulations are all done with forward evolution 
alone as they require less computational effort. 

The Born expansion of the echo operator up to second order reads 

Ma(0 = 1 - iAI{t) - A^}{t) + 0(a3), (2.9) 

with ^ ^ 

I{t)= f dTV{T), /(f) = / dT rdT'y(T)y(T') (2.10) 

JO Jo Jo 

and V'(f) = Uo{tyVUo{t) being the coupling in the interaction picture. 
Using this expansion we calculate the purity of the central system, aver- 
aged over the coupling and the Hamiltonian of the environment. 

The reader must notice that at no point we used that the dimension of the 
central system is 2. In fact, we only required the locality of the operator 
Uq and the fact that Ma(0 ~ 1. Thus, all this reasoning (in particular 



eqs. 2.8 2.9 and 2.101 is equally valid for a completely general spectator 



configuration or even more general configurations to be introduced later. 



2.3 The solution 

In appendix [A| we compute the average purity {P{t)) as a function of 
time in the linear response approximation eq. ( |2.9| , following the steps 
outlined in sec. 2.2 The average is taken with respect to the coupling Vi e 
[using eqs. E.4 and |E.5| , the random initial state \ipe), and the spectrum 
of He. In the limit of Ng oo, we obtain [eq. ( |A.8| ), eq. ( |A.36| )] 



(Pit)) 



1-2A^ dT dr' Re Aji(t,t') + 0(A4), (2.11) 



with 



Aji(t,t') = [Ci(|t-t'|)-Si(t-t')]C(|t-t'|)+^goe[1-SU-t- 



(2.12) 
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where Xgoe = 1 for the TRI case, and Xgob = for the non-TRI case. 
The correlation functions Ci{t),Si{t),S[{t), and C(t) are defined in ap- 
pendix A.5[ The first three depend on the state of the central system. 
Note that S[ (t) is only relevant in the case of a GOE, and curiously is not 
strictly a correlation function as it contains a dependence on the sum of 
both times. The last one, C(t) deserves special attention, since it depends 
on the spectral properties of the environment determined by the function 



1 



-iE jt 



C{t) = l + d{t/TH) 



(2.13) 



(recall eq. ( E.8| and subsequent equations). Here the Ey's are the eigenen- 
ergies of He and Th is the corresponding Heisenberg time. Actually the 
validity of eq. ( 2.11| is not dependent on the environment being repre- 
sented by a GOE (/6 = 1) or GUE (/3 = 2). For these the two-point form 

factor is well known l)Meh9H but any ensemble with the correspond- 
ing invariance properties will do, for example the FOE or FUE ||D^9TI. 

We first study the GUE case with and without an internal Hamiltonian 
governing the qubit. The next step is to work out the GOE case. There 
we concentrate on the case with no internal Hamiltonian governing in the 
qubit since we want to keep the discussion as simple as possible to focus 
on the consequences of the weaker invariance properties of the ensemble. 



2.4 The GUE case 



We are now in the position to give an explicit formula for in the 

GUE case. This formula will generally depend on some properties of 
the initial condition We wish to write it in the most general way. 

However the symmetries involved in the problem reduce the number of 
parameters needed to describe the initial state | t/^i ) . 

Recall that H = Hg + Hi + XV represents an ensemble of Hamiltonians in 
which He and V are chosen from GUEs of dimension Ne and 2Ne respec- 
tively, whereas Hi together with the initial condition \-ipi) remain fixed 
throughout the calculation. The operations under which the ensemble is 
invariant are local (with respect to the partitioning of the Hilbert space 
into Til and T^e)/ unitary (due to the invariance properties of the GUE), 
and leave Hi invariant. Hence the transformation matrices must be of the 



20 



Chapter 2. One qubit decoherence 



form 

LJOexp(iaHi) (2.14) 

with a a real number and U a unitary operator acting on T^e- The solution 
must also be invariant under that transformation. 

This freedom allows to choose a convenient basis to solve the problem. 
On the one hand, it allows to write Hq in diagonal form (as done during 
the discussion of appendix |A]), and on the other hand, we can use it to 
represent the initial state of the qubit in such a way that there is no phase 
shift between the two components of the qubit. This can be achieved by 
appropriately choosing a in eq. ( |2.14 |. We thus write, without loosing 
generality 

|t/;i) = cos(/)|0) + sin</)|l) , (2.15) 

where |0) and |1) are eigenstates of Hi. Notice that \i (p G {0, 7r/2}, 
is an eigenstate of Hi. Finally, we choose the origin of the energy 
scale in such a way that the Hamiltonian of the qubit can be written as 

Hi = (A/2)|0)(0|-(A/2)|l)(l|. 



We obtain the average purity from the general expression in eq. (2.11 1 and 
eq. ( 2.12| . For a pure initial state pi = \ipi){^i\ the relevant correlation 



functions ReCi(T), Si(t), and C(t) are given in eq. (A.46I, eq. (A.48I, and 



eq. ( A.41| |, respectively. Using the symmetry of the resulting integrand 



with respect to the exchange of t and r', we find 

(P(f)) = l-4A2|'dT At' C{t') [l-g^ (1- cos At')] +0{\\^-^) 

(2.16) 

with 

A .A 3 + cos(4(*) 
= cos^ cp + sin^ (p = ^^-^ (2.17) 

quantifying the "distance" between \^p) and an eigenbasis of Hi. 

Let us consider the following two limits for Hi. The "degenerate limit", 
where the level splitting A is much smaller than the mean level spac- 
ing de = 2n/ Th of the environmental Hamiltonian, and the "fast limit", 
where the level splitting is much larger. In the latter case, the internal 
evolution of the qubit is fast compared with the evolution in the environ- 
ment. (We shall refer to these limits also in later sections.) 

The degenerate limit leads to the known formula [|PS07| 

Pd(0 = 1-AVth(0. (2-18) 
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with 

fr^it) = 2fmax{f,TH} + -^(min{i,TH})^. (2.19) 

The result does not depend on the initial state of the qubit. Due to the 
degeneracy all states are eigenstates of Hi and thus equivalent. The lead- 
ing term of the purity decay is linear before the Heisenberg time and 
quadratic after the Heisenberg time. Similar features were already ob- 
served in fidelity decay and purity decay in other contexts |iGPSZ06l . 



In the fast limit (A » dg), purity is obtained from eq. (2.161 by replacing 
cos At' by 1 when it is multiplied with the S function [see eq. ( 2.13| ], and 
by zero everywhere else. For finite Ng care must be taken, since we are 
assuming Zeno time (which is given by the "width of the (5-f unction") to 
be much smaller than all other time scales, such that A <C Ned^. The 
resulting expression is 

Pf(0 = 1 - A2[(1 -^^)/.h(0 +2^<?.^Th] . (2.20) 

Typically (unless \tpi) is an eigenstate), this formula displays a domi- 
nantly linear decay below the Heisenberg time, and a dominantly quadratic 
decay above, similar to eq. ( 2.18| . 



In fig. |2.l| we compare numerical simulations of the average purity {P{t)) 
(symbols) with the corresponding linear response result (dashed lines) 
based on eqs. eq. ( 2.18| and eq. ( 2.20[ |. The numerical results are obtained 



from Monte Carlo simulations with 15 different Hamiltonians and 15 dif- 
ferent initial conditions for each Hamiltonian. We wish to underline two 
aspects. First, the energy splitting in general leads to an attenuation of 
purity decay. Even though a strict inequality only holds for the limiting 
cases, P^it) > Puit) (for t ^ 0), we may still say that increasing A tends to 
slow down purity decay. This result is in agreement with earlier findings 
on the stability of quantum dynamics [PZ0 2J. Second, for the fast limit 
and an eigenstate of Hi (g^ = 1) we find linear decay even beyond the 
Heisenberg time. A similar behaviour has been obtained in ||GPSZ06t , but 
there an eigenstate of the whole Hamiltonian was required. 

In IIPS07I it was shown that exponentiation of the linear response result 
leads to very good agreement beyond the validity of the original approx- 
imation. We use the formula eq. ( |C.1| 



Pelr(0 =-P~ + (l--Poo)exp 



1 - PLRjt) 
1-Pco 



(2.21) 
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Figure 2.1: Numerical simulations for the average purity of one qubit 
using non-TRI Hamiltonians, as a function of time in units of 
the Heisenberg time th of the environment. For the coupling 
strength A = 0.01 and Ng = 2048, we show the dependence of 
{P{t)) on A (the level splitting) and (p (characterizing the ini- 
tial state) of the internal qubit: A = (green circles), A = 8, 
(p = 7z/A (blue squares), and A = 8, (p = (red triangles). 
The corresponding linear response results (dashed lines) and 
exponentiated linear response results (solid lines) are based on 
eq. ( 2.16| and eq. (2.211, where Poo is given in|C] Note that the 
level splitting A tends to slow down decoherence. 
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where PlrIO is truncation to second order in A of the expansion eq. (2.16 i, 
and Poo = 1/2 the estimated asymptotic value of purity for t ^ oo, see 
appendix|C] From fig. 2.1 we see that the exponentiation indeed increases 
the accuracy of the bare linear response approximation. 



2.5 The GOE case 



We drop Hi for a moment, leaving Hq = He, resulting in 

Ha = He + AV. (2.22) 

He is chosen from a GOE of dimension Ng and acts on Ti.^,; V is chosen 
from a GOE of dimension 2Ne and acts on 'He ^ ^i- The resulting ensem- 
ble of Hamiltonians is invariant under local orthogonal transformations. 
In the qubit this invariance allows rotations of the kind exp(iat7'i^) G 0{2). 
If such transformations are represented on the Bloch sphere, they become 
rotations around the y axis. Hence, they can take any point on the Bloch 
sphere onto the xy-plane. Supposing this point represents the initial state, 
it shows that we may assume the initial state in the qubit to be of the form 

I*,) ^ (2.23, 

In this expression, 7 G [—7z/2, re/ 2] denotes the angle of the vector rep- 
resenting the initial state with the ;cz-plane (see fig. |2.2||. 



In order to obtain the linear response expression for {P{t)) we again make 
use of eq. ( 2.11| and eq. ( 2.12[ |. However, apart from the correlation func- 



tions used in the GUE case, we have now to consider in addition S'^(t), 
as given in eq. ( A.53| . The special case Hi = can be simply obtained by 



setting A = 0. This yields 

Aji(t,t') = C(|T-T'|)+sin2 7. (2.24) 
After evaluating the double integral in eq. ( 2.11| |, we obtain 



where 



{P{t)) = l-A^{e[3-cos{2j)]+2tTH-2B'^\t)Y (2.25) 
B^^Ht) =2^*dT^^dT'4^^(T7TH) (2.26) 
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Figure 2.2: Any pure initial state of the qubit can be mapped onto the 
Bloch sphere. Here, we show the angle 7 defined in eq. ( 2.23| in 
color code. Regions of a given color represent subspaces which 
are invariant under the transformation exp(iac7'y). 



is the double integral of the form factor. It can be computed analytically, 
but the resulting expression is very involved IIGPS04| . For our purpose 

it is sufficient to note that for t <C Th, B^\t) a (as in the GUE case), 

whereas for t 3> Th, t — B2\t) grows only logarithmically. Thus it has a 
similar behavior as in the GUE case, eq. ( G.6| . 



In fig. 2.3 we show {P{t)) for 7 = (green squares), for 7 = Ti/2 (blue 
circles), and for random values in the whole Bloch sphere (red triangles). 
In contrast to the GUE case in the degenerate limit, the average purity 
depends on the initial state (via the angle 7). The fastest decay of purity is 
observed for 7 = n/2, where the image under the time reversal operation 
becomes orthogonal to the initial state. The slowest decay is observed for 
7 = 0, which characterizes states which remain unchanged under the 
time reversal symmetry operation. In the lower panel we show numerical 
results for the standard deviation of the purity as a function of Ng, the 
dimension of the Hilbert space of the environment. We consider the same 
cases as on the upper panel: random initial states \ipi) with fixed 7 = 
(green squares), with fixed 7 = 7r/2 (blue circles) and random states 
\xp-[) uniformly distributed on the Bloch sphere (red triangles). Note that 
along with the random choice of \tpi), also Hg, V-[ g, and \xpe) are randomly 
chosen from their respective ensembles. We clearly see that for those 



2.5. The GOE case 



25 





Figure 2.3: We display the behavior of purity and von Neumann en- 
tropy, in different regions in the Bloch sphere connected by or- 
thogonal transformations, characterized by 7 in eq. ( 2.23| . On 
the top figure we see the enveloping curve after running 100 ini- 
tial conditions (thick regions), their average (symbols) and the 
predicted behavior (solid curves) by eq. (2.251. If 7 = we use 
color green; if 7 = 7r/2 we use blue; and if we allow arbitrary 
7 we use red. In the lower panel a is plotted for t = 40. We use 
the same coding as the upper figure. For a fixed value of 7 (blue 
and green) there is asymptotic self averaging whereas for an ar- 
bitrary initial condition (red) the standard deviation reaches the 
finite value predicted in eq. ( 2.28|, p lotted as a horizontal line. 
We fixed A = lO^"'. A line oc l/i/Ne is also included. 
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-1/2 
e 



cases where 7 is kept fixed, the standard deviation falls off like N, 
By contrast, the standard deviation converges to a finite value, when | tpi ) 
chosen with no restriction. That value can be estimated from the standard 
deviation of the function cos 27, which yields: 



Tz/2 



nil L 



cos(27) — cos(27) 



2 cos(27) 



d7 



16 
45' 



(2.27) 



Assuming that for Ne — > 00, the fluctuations of cos 27, is the only source 
for the fluctuations of purity, the standard deviation of the purity is 



3^5 



4 Ar-l^ 
e I 



(2.28) 



The statements of the last two paragraphs can be directly translated to 
the von Neumann entropy S. For one qubit it has a one to one relation 
with purity 

S{P) = h I ) + ^ I 2 ) ' " -x\og2X. 

(2.29) 

Observe the entropy scale on the upper panel of fig. 2.3 The conse- 
quences of the weaker invariance properties of the GOE, and the relation 
to the states will be analyzed in detail in a general framework in a later 
paper. 



Chapter 3 

Two qubit decoherence 



In this chapter, we address the question whether entanglement within 
a given system affects its decoherence rate. In particular, as the name 
suggests, we are going to study two qubit decoherence. We will use 



the spectator model described in sec. 1.2 Moreover we shall consider 
the first nontrivial example thereof, in the sense that the spectator space 
plays a nontrivial roll. However it is still the simplest example allowing 
this possibility as both the spectator and the interacting space are qubits. 
We shall study two other configurations, namely when both qubits are 
coupled to one or two environments. There, we shall start appreciating 
the power of the spectator model; we are going to be able to express there 
decoherence in terms of the decoherence in the spectator configuration. 

We will base our arguments on the calculations made in appendix [A| and 
the results discussed in the previous chapter. Again we are going to 
study the linear response regime, and test with Monte Carlo simulations 
the validity of a heuristic exponentiation. The symmetries of the clas- 
sical ensembles will continue to play an important role to simplify the 
problem. Moreover we shall find that entanglement has the property of 
transporting the symmetry from one qubit to the other. 



We first describe the models (configurations) we are going to study, sec. 3.1 



Next we analyze in detail the results for the simplest configuration, the 



spectator model for two qubits, sec. 3.2 There we consider both the GUE 



and GOE cases. We then generalize the result to the other configurations 
in sec. 13.31 and sec. 13.41 
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Figure 3.1: Schematic representations of the different dynamical con- 
figurations studied in this article: (a) the 2 qubit spectator 
Hamiltonian, (b) the separate, and (c) the joint environment 
Hamiltonian. 



3.1 The models 



For the two qubit case, the Hilbert space structure must be slightly more 
complicated than eq. (2.1 1. We need at least to provide the Hilbert space 
for a second qubit, and, in one of the models, we shall need an additional 
Hilbert space for a second environment. In all cases the initial condition 
is pure in the central system, but the two qubits can share some entangle- 
ment. We shall consider three different dynamical scenarios, all explicitly 
excluding any interaction between the two qubits: 



(a) The 2 qubit spectator Hamiltonian: 
The Hilbert space structure is 

n = ni0n20ne. (s.i) 

Both Hi and Hi are the state spaces of the qubits (dim Hi = dim H2 = 
2) and He (of dimension Ne) is the state space of an environment. 
The central system is, obviously, 

Hc = Hi® Hi. (3.2) 

Only the first qubit is coupled to an environment, and we allow for 
local d}mamics. The total Hamiltonian reads 

Ha = Hi + H2 + He + AVi^e , (3.3) 



where A is again a real number modulating the strength of the cou- 
pling. We recall the remark in sec. 1.2 namely that the sub-indices 
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of the operators indicate the space on which they act. This situation 
is shown schematically in fig. |3.1[ a). Notice that if we choose an ini- 
tial state where the two qubits are already entangled, this provides 
the simplest situation which allows to study entanglement decay. If 
the two qubits are initially not entangled, the process reduces ef- 
fectively to the one-qubit case, described previously sec. 2.1 The 
special case Hi = H2 = has been considered in Ref. IIPS07L 

The reader must notice that this situation is an example of the spec- 



tator configuration, discussed in sec. 1.2 In this case the interacting 



subspace is Hi and the spectator is played by 7i2- 

(b) The separate environment Hamiltonian: 
The Hilbert space structure is 



n = ni®n2®ne® n^'. (3.4) 



Besides the subspaces in eq. ( |3.1| (which keep their meaning), we 
have an additional space "Hg' which represents a new environment. 
We again allow for similar dynamics except that we allow the sec- 
ond qubit to interact with the new environment. The two environ- 
ments are assumed to be non-interacting and uncorrelated: 

\^{^)) = \^ll)\^e)\%'), eUi^Hi, \^e) G7i:e, and|l/^e') ^ H,, 

(3.5) 

Thus, the Hamiltonian of this model reads 

Ha„a2 = Hi + H2 + He + He. + Ai V^^, + A2 V^,^, , (3.6) 

where V2,e' and He' describe the coupling to - and the dynamics 
in the additional environment. Both quantities are chosen inde- 
pendently from the respective random matrix ensembles, in perfect 
analogy with Vi e and H^. The real parameters Ai and A2 fix the cou- 
pling strengths to either envirorunent. This model, see fig. |3.1| ^b), 
may describe two qubits that are ready to perform a distant telepor- 
tation, where each of them is interacting only with its immediate 
surroundings. It can also represent a pair of qubits that, although 
close to each other, interact with different and independent degrees 
of freedom. 

(c) The joint environment Hamiltonian: 



The third case, shown in fig. 3.1 'c), describes a situation in which 



both qubits are coupled to the same environment, even though the 
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coupling matrices are still independent. The Hilbert space is iden- 
tical to the one for the 2 qubit spectator corvfiguration. The total 
Hamiltonian reads 

Hx,M = Hi + H2 + He + Ai yi,e + A2 V2,e , (3.7) 

where V2,e describes the coupling of the second qubit to the envi- 
ronment. It is chosen independently from the same random matrix 
ensemble as Vi^e- 



3.2 The spectator Hamiltonian 



The first step to calculate the decoherence of the initial state 

^0 = |'/'12)(V'12| O \^pe){^e\, 



(3.8) 



evolved with the Hamiltonian ( |3.3| , is to realize that the echo operator 
does not contain H2. The quantum echo of ^0 after time t is 



>Ml{t) 



Since Q (t) remains a pure state in T-Li®'H2® 'He, 

P{t) =trpc{tf = trpe{tf 



(3.9) 



(3.10) 



with pc{t) = tVeQ^it) and pe{t) = tVcQ^it)- This simply means that, 
as a formality, we can calculate purity of the central system, calculating 
purity of the environment. As the echo operator acts as the identity on 
the second qubit, 

pe{t) = triMA(0(tr2(?o)Ml(0 (3.11) 

= triMA(f) ipi(^\lpe){^e\)Ml{t), 
where pi = tr2 \ tpi2){^i2\- 

We may therefore compute the purity of the spectator model, without 
ever referring explicitly to the second qubit! Any dependence of the decay 
of purity on the central system as a whole is encoded into the initial 
density matrix pi . This also implies that we can use the results obtained 
in [a} and hence eq. ( 2.11| and eq. ( 2.12| remain valid. The only difference 
is that for the correlation functions Ci(t), Si(t), and S[{t), we now have 
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to insert the respective expressions which apply for mixed initial states 
of the first qubit. These expressions are given in |A.5 We stress, for later 
reference, that this line of reasoning is not limited by the fact that the 
interacting and spectator systems are qubits. 



3.2.1 The GUE case 

We wish to write the initial condition in its simplest form. We must 
respect the structure of the Hamiltonian ( |3.3| |. However we can still take 
advantage of all its invariance properties, when seen as an ensemble. 
Given a fixed Hj, that ensemble of Hamiltonians is invariant under local 
operations of the form 

® exp iocHi (g) U2 (3.12) 

where U^^ G U{Ne) is any unitary operator acting on the environment, 
OL a real number, and U2 G U{2) is any unitary operator acting on the 
second qubit. 

The freedom within the qubits allows to choose a basis {|0),|1)}(8){|0),|1)} 



in which the initial state [see eq. (3.8 1] can be written as 



\rpi2) = cos0(cos^|O) +sin^|l))|0) +sin0(sin(^|O) -cos(/)|l))|l), (3.13) 

and still. Hi = ^ |0) (0| — f |1) (1 1 is diagonal. Let us prove it in a construc- 
tive manner. To find this basis we start using the Schmidt decomposition 
to write 

\1p12) = COS0IO1O2) +sin0|iii2) (3.14) 

with {|0j), |1,)} being an orthonormal basis of particle i. For the first 
qubit, we fix the z axis of the Bloch sphere (containing both |0) and |1)) 
parallel to the eigenvectors of Hi, and the y axis perpendicular (in the 
Bloch sphere) to both the z axis and |0i). The states contained in the xz 
plane are then real superpositions of |0) and |1), which implies that 

|0i) = cos^lO) + sin^|l), (3.15) 
= sin^lO) -cos^|l), (3.16) 

for some (p. In the second qubit it is enough to set 

IO2) = |0), (3.17) 
II2) = |1). (3.18) 
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Figure 3.2: A figure to visualize the way the initial condition is 
parametrized, using the Bloch sphere, is presented. On the left, 
qubit one has an internal Hamiltonian. Its eigenvectors (|0) and 
|1)) are represented in blue. The z axis is chosen parallel to |0). 
The X axis is chosen so that both |0i) and have real coeffi- 
cients; i.e. such that the xz plane contains |0i) and On the 
right we represent the second qubit. We have absolute freedom 
to choose the basis (even if an internal Hamiltonian is present), 
and thus we choose it according to the natural Schmidt decom- 
position. 



This freedom is also related to the fact that purity only depends on 



tr2 IV^ii) (V'lil- A visualization of this argumentation is found in fig. 3.2 
The angle 9 G [0, 7r/4] measures the entanglement 

C(|i/'i2)(t/'i2|) =sin20 (3.19) 

whereas the angle ^ G [0, 7r/2] is related to an initial magnetization. 
The general solution for purity using this parametrization is 

P{t) = 1 - 4A2 dT dT'C(T') + gfl COS At'] + 0(A^ N'^), 

(3.20) 

where the geometric factors gf^^ G [0,1/2], and g^^^ G [1/2,1] are ex- 
pressed as 

S^el =ge{l-g^)+g^{l-ge) (3.21) 
gfl =2{l-ge)-g^{l-1ge), (3.22) 
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Figure 3.3: Visualization of the geometric factors g^^^^, and g^"^ from 

left to right respectively. For ^^^^ we see that for pure eigenstates 
of Hi its value is zero. This leads to a higher qualitative stability 
of this kind of states. 



in terms of the functions g(p and ge, defined in eq. (2.171. Both geometric 
factors are shown in fig. |3.3[ Eq. ( 3.20| is obtained from eq. (2.111 and 
eq. ( 2.12| by insertion of the eq. (A.46I and eq. ( A.52| for ReCi(T) and 
Si(t), respectively. 

We consider again two limits for A. In the degenerate limit (A ^ 1/th) 
purity decay is given by 



■ge) fruit)' 



(3.23) 



where /^h (t) is defined in eq. ( 2.19| . The result is independent of (p since a 
degenerate Hamiltonian is, in this context, equivalent to no Hamiltonian 
at all. The 0-dependence in this formula shows that an entangled qubit 
pair is more susceptible to decoherence than a separable one. 

In the fast limit (A » 1/th) we get 



Pf(0 = 1 - 



(;^/.„(o+2TH^gf 



(3.24) 



For initial states chosen as eigenstates of Hi we find linear decay of purity 
both below and above Heisenberg time. In order for pi to be an eigenstate 
of Hi it must, first of all, be a pure state (in 'Hi). Therefore this behavior 
can only occur if 6 = or 9 = tt/2. Apart from that particular case, 
we observe in both limits, the fast as well as the degenerate limit, the 
characteristic linear/ quadratic behavior before/ after the Heisenberg time 
similar to the one qubit case. 
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Figure 3.4: Numerical simulations for the average purity as a func- 
tion of time in units of the Heisenberg time of the environment 
(spectator configuration, GUE case). For the coupling strength 
A = 0.03 we show the dependence of {P{t)) on the level split- 
ting A in Hi and on the initial degree of entanglement between 
the two qubits (in all cases (p = 7r/4). For separable state {9 = 0) 
and a degenerate limit we use whereas for the fast limit, with 
A = 8, we use A; For Bell states, the corresponding limits are 
encoded as • in the degenerate limit and ♦ in the fast limit 
(A = 0.8). The corresponding linear response results (dashed 
lines) are based on eq. ( |3.23 l, eq. ( 3.24| and eq. ( 2.21[ l. The expo- 
nentiated linear response results (solid lines) are obtained based 
on the results of appendix |C] The theoretical curves are plotted 
with the same color, as the respective numerical data. In all 
cases Ne = 1024. 



In fig. 3.4 we show numerical simulations for {P[t)). We average over 
30 different Hamiltonians each probed with 45 different initial condi- 
tions. We contrast Bell states (cp = n/A, 6 = n/ 4) with separable states 
{(p = n/A, 6 = 0), and also systems with a large level splitting (Ne 3> 
A = 8 ^ 1/th) with systems having a degenerate Hamiltonian (A = 0). 
The results presented in this figure show that entanglement generally en- 
hances decoherence. This can be anticipated from fig. 3.3 since for fixed 

.(1) 



(p, increasing the value of 6 (and hence entanglement) increases both 

(2) 

and gg^. At the same time we find again that increasing A tends to re- 
duce the rate of decoherence. A strict inequality only holds among the 
two limiting cases (just as in the one qubit case): From gf'^ = 2 — gf^^ — ge, 

it follows that (Pp - Pd)/A^ = gfl[fTn{i) - 2^Th] > 0. Therefore, for fixed 



3.2. The spectator Hamiltonian 



35 



initial conditions and f greater than 0, Pf(0 > Pvit)- This is the second 



aspect illustrated in fig. 3.4 



In order to extend the formulae to longer times/ smaller purities we expo- 
nentiate them using the results of appendix |C] The numerical simulations 
(see fig. |3.4| agree very well with that heuristic exponentiated linear re- 
sponse formula. In one case (blue rhombus) where the agreement is not 
as good, we found that it is the inaccurate estimate Poo = 1/4 of the 
asymptotic value of purity, which leads to the deviations. 



3.2.2 The GOE case 



Let us consider the GOE average of both He and Vi g. When averaging 
He and V-[ g over the GOE, we are again confronted with the fact that the 
invariance group is considerably smaller than in the GUE case. In this 
context the initial entanglement between the two qubits has a crucial im- 
portance since it "transports" the invariance properties from the spectator 
to the coupled qubit. 

For the sake of simplicity, we focus on the degenerate limit setting Hi = 0. 
Note that on the basis of the results in appendix |A| the general case can 
be treated similarly and the corresponding result will be presented at the 
end of this subsection. 

We first specify the operations under which the spectator Hamiltonian 



eq. (3.3 1, considered as a random matrix ensemble, is invariant. As both, 
the internal Hamiltonian of the environment and the coupling, are se- 
lected from the GOE the invariance operations form the group 

0(Ne) X 0{2) X U{2), (3.25) 

and have the structure O^^ exp{mcry) (g) U2, with O^^ being an orthog- 
onal matrix (acting in He), a a real number, and U2 a unitary operator 
acting on the spectator qubit. The direct product structure of the invari- 
ance group obliges us to respect the identity of each particle, but allows 
to analyze each qubit separately. For instance, if we replace the random 
coupling matrix Vi e with one which involves both qubits, the invariance 
group would be 0(Ne) x 0(4). As a consequence, purity decay would 
become independent of the entanglement within the qubit pair: for any 
entangled state one can find an orthogonal matrix which maps the state 
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onto a separable one 



We write the initial condition \1p12) as in eq. (3.14|. For the coupled 
particle follows the same analysis made in sec. 2.5 We can thus write 
|0i) = 2^^^-^[\0) +exp(i7)|l)] and, in order to respect orthogonality, 
= 2^^''^exp(i^)[|0) — exp(i7)| l)]. F or the second qubit we have the 
same complete freedom as in eq. ( 3.2. We thus select IO2) = |0) and 
II2) = exp(— i^) |1) to erase the relative phase in the first qubit and finally 
write the initial state as 



cos0(|O) +£'iT|l))|0) +sin0(|O) -e'T|l))|l) 



(3.26) 



The average purity is still given by the double integral expression in 
eq. ( 2.11| . However, in the present case the mixed initial state pi = 
tr2 1 1/^12) {^12] must be used. For A = 0, the resulting integrand reads 

Aji(t, t') = (2 - ge)C{\T -r'\) + l-ge + {2ge - 1) sin^ 7, (3.27) 



where C ( | t 
we obtain 



t'I) is given in eq. (2.131. Evaluating the double integral. 



2cos2(20)cos2 7] + (4-2^0) 



,(1)/ 



where B^^^(f) is given in eq. (^2. 26b 



(3.28) 

As in the GUE case, this result de- 
pends on the entanglement of the initial state and, as in the one-qubit 
GOE case, it also depends on 7. Again it turns out that Bell states are 
more susceptible to decoherence than separable states. Note however, 
that purity as a function of d is not monotonous. Hence, a finite increase 
of entanglement does not guarantee that the purity decreases everywhere 
in time. For separable states, gQ = 1, we retrieve formula eq. (2.251. How- 
ever, for completely entangled states, d = n/A the dependence on 7 is 
lost. This is understood from a physical point of view, noticing that 
any local unitary operation on a Bell state can be reduced to a single 

-^One can see a general state \^ri) = 12^=o{'^i + i^jOlO' "j'^j € K as characterized by 
2 real vectors a = {uq, a^, 02,(13,), b = 0{bQ,bi,b2,b3) 6 R*. We can rotate « to a vector 
having only the first component a' = (flg, 0,0,0) using an orthogonal matrix O. Consider 
the vector b' = Ob rotated by the same transformation. We next repeat the procedure on 
the last 3 axes (for b') with an orthogonal transformation O' , to obtain a vector having 
only the first two components. We thus have, an orthogonal transformation that takes 
\xp12) to 0'0\xpi2) = Li=o{ai +ib,")\i) = |0i)(a|02) +^|l2), a separable state. 
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local unitary operation acting on a single qubit. hi other words, given 
Ui^2 = Ui® U2, there exists a unitary U[ such that 

Ui O LJ2|Bell) = U[ ® l|Bell) (3.29) 

(|Bell) is any 2 qubit pure state with C = 1, e.g. |00) + We can then 
say that the invariance properties in the coupled qubit are inherited from 
the spectator qubit via entanglement. 

Let us obtain the standard deviation for the different possible initial con- 
ditions in the qubits. We want to analyze the situation separately for a 
fixed value of concurrence. Then, as the invariant measure of the ensem- 
ble of initial conditions, and fixing the amount of entanglement, we shall 
use the tensor product of the invariant measures in each of the qubits. 
Since there is no dependence of eq. ( 3.28| on the second qubit, the appro- 



priate invariant measure is trivially inherited from the invariant measure 
for a single qubit. The resulting value for the standard deviation is 

= 4M!!f5^. (3.30) 
Sv 5 

Based on the appendix |A] we can also obtain the average purity for A / 
0. The parametrization of the initial states is more complicated since 
two preferred directions arise, one from the eigenvectors of the internal 
Hamiltonian and the other from the invariance group. The result can be 
expressed in the form given in eq. ( 2.11| , with 



ReAji(T,T') = C(|T-T'|) [^g + ^g cosA(t-t')] 

+ 8e]l - (1 - ^S) cos[A(T + t') - 2^] + 0{A\ N-') . 

The angle ?/ is related to a phase shift between the components of any of 
the eigenvectors of the initial density matrix pi. Here we can again see 
the term containing a sum of times, i.e. a term that is not a correlation 
function. However, this term is small, and adequate values must me 
chosen to see its effect. Still it is observable in numerical simulations, 
with moderate effort. 



3.3 The separate environment Hamiltonian 

We proceed to study purity decay with other configurations of the envi- 
ronment. Consider the separate environment configuration, pictured in 
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fig. 3.1 h). The corresponding uncoupled Hamiltonian is 

Ho = Hi + H2 + He + Hgf 

and the coupHng is 



(3.31) 



AV = AiVe,! + A2Ve',2- (3-32) 

From now on we assume that the internal Hamiltonians of the environ- 
ment and the couplings are chosen from the GUE. The generalization 
for the GOE can be obtained along the same lines using the correspond- 
ing results of sec. 3.2.2 The initial condition has a separable structure 
with respect to both environments, see eq. ( |3.5| . The coupling in the in- 
teraction picture separates into two parts acting on different subspaces 
Ay = Aiy(i) + \2V^^\ where 

t/(l) = ei(Hi+He)y^^^g-i(Hi+He)^ y(2) ^ ^i(H,+H^,) y^^ ^^^-i{H,+H^,) _ (333) 

Notice that y'^) (t/^^)) does not depend on H^, (He). Since y(^' and y'^) 
are uncorrected, quadratic averages separate as 



K2)r/(2)N 



(3.34) 



This leads to a natural separation of each of the contributions to purity 



1 - (p(0) = 1 - p^iUt) + 1 - p^pUt), 



,(2) 



(3.35) 



where Pspec(0 denotes the average purity with particle i being a spectator, 
as given in sec. 3.2 In this way, the problem formally reduces to that of 



the spectator model. The respective expressions in sec. 3.2 may be used. 
For instance, if we assume broken TRI, we obtain from eq. (3.201 



(2) 



+ ^ri C0SA;T' Q(t') +0(A4,Ne-l), 

(3.36) 

where Ci and C2 are the correlation functions of the corresponding envi- 
ronments defined in exact correspondence with eq. ( 2.13| l, for Hg and He' 
respectively. If in one or both of the qubits, the level splitting in the inter- 
nal Hamiltonians is very large /small compared to the Heisenberg time 
in the corresponding environment (denoted by Te and Te' for He and Hgv 
respectively) the degenerate and /or fast approximations maybe used. As 
an example, if Ai <C 1 / Te and A2 <C 1 / Tg' we find 



Pvit) 



{2-g8){^lfrM+4fr,m 



(3.37) 
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whereas if Ai » 1/Te and A2 » 1/^ 

(3.38) 

It is interesting to note that if we have two separate but equivalent envi- 
ronments (i.e. both Heisenberg times are equal), we get exactly the same 
result as for a single environment. Also notice that the Hamiltonian of 
the entire system separates and thus the total entanglement of the two 
subsystems {Hi ^ He and H2 H^') becomes time independent. 



Pvit) 



A? 



t) + 2Te 



4 



3.4 The joint environment configuration 



The last configuration we shall consider is the one of joint environment; 

(3.39) 



see fig. 3.1 'c). Its uncoupled Hamiltonian is 

Ho = Hj + H2 + He 
whereas the coupling is given by 



(3.40) 



Notice the similarity with eq. (3.31 1 and eq. (3.32 1. However, as discussed 



in the introduction, they represent very different physical situations. The 
coupling in the interaction picture can again be split AV = A^V^^^ + 
A2V^^\ where 

^ gi(Hi+He)y^^e-i(Hl+He)^ y(2) ^ gi(H2+He)y^ ^g-i(H2+He) _ (34;^) 



Note the slight difference with eq. (3.33 1. However y(^) and V^^^ are still 
uncorrected, enabling us to write again eq. ( 3.34| |. 



From now on, the calculation is formally the same as in the separate 
environment case. Hence we can inherit the result eq. ( 3.36| directly, 
taking into account that since they come from the same environmental 
Hamiltonian, the two correlation functions are the same. In case any 
of the Hamiltonians fulfills the fast or degenerate limit conditions, the 
corresponding expressions to eq. ( 3.37| and eq. ( 3.38| can be written. As 
an example, if the first qubit has no internal Hamiltonian, and the second 
one has a big energy difference, the resulting expression for purity decay 
is 
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(3.42) 
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where th is the Heisenberg time of the joint environment. Monte Carlo 
simulations showing the validity of the result were done with satisfactory 
results, comparable to those obtained in fig. 3.4 The parameter range 
checked was similar to that in the figure. 



Chapter 4 

Entanglement decay 



In the previous chapter we studied purity decay of two qubits. Purity 
measures entanglement with the environment, but we can wonder how 
decoherence affects the internal quantum mechanical properties of the 
central system. Possibly the most important quantum mechanical prop- 
erty of a multi-particle system is its internal entanglement. As discussed 
in appendix |D| a simple and meaningful measure of two qubit entangle- 
ment is concurrence (C). Since concurrence is defined in terms of the 
eigenvalues of a Hermitian 4 x 4-matrix, an analytical treatment, even in 
linear response approximation, is much more involved than in the case of 
purity. A study along the same lines followed in the last chapters, is out 
of reach for the time being. 

We shall first explore a relation (first found in fPS06], partly explained in 
IIZBOS I, and further studied in [PS0 7I) between concurrence and purity. 
We show that this relation is valid in a wide parameter range (sec. |4.1| . 
Combining it with an appropriate formula for purity decay, we obtain an 



analytic prediction for concurrence decay in sec. 4.2 We compare our pre- 
diction with Monte Carlo simulations. It is essential that the two qubits 
do not interact; otherwise the coupling between the qubits would act as 
an additional sink (or source) for internal entanglement - a complication 
we wish to avoid. In the last part of the section we extend our ideas to 
non-Bell states. 
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4.1 The CP plane 



We study the relation between concurrence and purity using the CP- 
plane, where we plot concurrence against purity with time as a parame- 
ter. This plane is plotted in fig. 4.1 The gray area indicates the region of 



physically admissible states [ jMJWKOll . The upper boundary of this re- 
gion is given by the maximally entangled mixed states. This set depends 
on the entanglement measures chosen |WNG"''03| . Large purity i.e. low 
entanglement with the environment, is required in order to have large 
concurrence, i.e. entanglement within the pair. This is a consequence 
of the monogamy of entanglement. Note also that concurrence becomes 
identically zero before purity reaches its minimum at 1/4. Another re- 



gion of interest, plotted in red in fig. 4.1 'a), corresponds to those states 
(density matrices), which form the image of a Bell state under the set of 
local (i.e. acting separately on each qubit) unital operations liZBOS] . Uni- 
tal operations are those which preserve the identity |Key02( . Single qubit 
unital operations include bit flip and phase flip, whereas an example of a 
non-unital operation is amplitude damping l|NCOOL Finally, the Werner 

states pw = + (1 — a) |Bell) (Bell|, < a < 1, define a smooth curve on 
the CP-plane (black solid line). The analytic form of this curve is IIPS07I 

Cw(P) = max |o Vl2P-3-l |^ ^^^^ 

and will be referred to as the Werner curve. Note that states mapped to 
the Werner curve are not necessarily Werner states. 



In fig. 4.1 'b), we perform numerical simulations in the spectator configu- 
ration assuming broken time reversal symmetry (GUE case). We compute 
the average concurrence for a given interval of purity using 15 different 
Hamiltonians and 15 different initial states for each Hamiltonian. We fix 
the level splitting in the coupled qubit to A = 1 and consider two different 



values A = 0.02 and 0.14 for the coupling to the environment. Fig. 4.1 'b) 
shows the resulting CP-curves for different dimensions of Tie- Observe 
that for both values of A, the curves converge to a certain limiting curve 
as dim('He) tends to infinity. While for A = 0.02, this curve is at a finite 
distance of the Werner curve, for A = 0.14 it practically coincides with 
Cw(P)- Varying the configuration, the coupling strength, the level split- 
ting, or the ensemble (GOE/GUE), gives the same qualitative results in 
the CP plane, for large dimensions. In some cases we have an accumu- 
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Figure 4.1: We present the CP plane. In (a) we show the area of 
concurrence-purity (CP) combinations which are allowed for 
physical states (gray area plus the set {{0,P),P G [1/4,1/3]}). 
The image of a Bell state under the set of local unital operations 
defines the red area. The Werner curve eq. (4.1 1 is shown as a 
thick black solid line. In (b) we show curves {{C{t)), {P{t))) as 
obtained from numerical simulations of the spectator Hamilto- 
nian. GUE matrices are used during all numerical experiments 
in this section. We choose A = 1 and vary dim('We) for two 
different coupling strengths A = 0.02 (thick lines) and A = 0.14 
(thin lines). The resulting curves are plotted with different col- 
ors, according to dim(7^e) as indicated in the figure legend. 
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lation towards the Werner curve, in others there is a small variation, but 
staying in the unital area. 

Are the channels induced by our RMT model (RMT channels) typically 
unital? If so, we could use the results of [ZB05 I to explain our numerical 
findings. The fact that our RMT channels map Bell states to the unital re- 
gion only is not enough to ensure unitality. Non-unital channels can also 
map Bell states into the unital region. The amplitude damping channel 
acting on a single qubit provides the simplest example. We now examine 
the question in more detail. 

For small dimensions the channels are clearly not unital as can be inferred 



from fig. 4.2 Bell states are mapped, after the action of our RMT chan- 
nels, to the non unital region. The probability of these events, however, 
decrease rapidly with the environment size. Within numerical accuracy, 
the probability of being mapped outside the unital region, for 5 = 0.2 
and a time range between and 4000 is ~ 10^^ already for Ng = 32. For 
larger dimensions a tighter test must be performed as virtually all points 
in the CP plane are in the unital region. 

We want to test how close to unitality the RMT channels iS'rmt are. The 
channels depend on the particular Hamiltonian chosen from the ensem- 
ble, on time, and on the initial condition in the environment. As we want 
to test unitality of the channel on a single qubit, we can think of the spec- 
tator configuration, or the single qubit model. Since we want check how 
closely £^rmt(1) = 1/ we prepare an initial pure condition (in the qubit 
plus environment) that leads to a completely mixed state in the qubit, i.e. 

|0#) + |1-/-P^) 

with {^g^\ipg ^) = Sjj. We let the state evolve with a particular member of 
the ensemble of Hamiltonians defined in eq. ( |2.2| . Afterwards we eval- 
uate the Euclidean distance d{-,-) in the Bloch sphere, of the points cor- 
responding to the resulting mixed state in the qubit and the fully mixed 
state. For unital channels this distance should be exactly zero. The av- 
erage distance is plotted as a function of the size of the environment in 



fig. 4.3 for two coupling values, and various times. Instead of reporting 
the times, we report the approximate value of concurrence a Bell pair 
would have after the corresponding time, and thus the area to which it 
would be mapped in the CP plane. We conclude form the figure that 
the unitality condition is approached algebraically fast as the size of the 
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Figure 4.2: We show the points reached in the CP plane after the evo- 
lution with different sizes of the environment. Four initial con- 
ditions (with different colors) are studied per environment size. 
For very small environments the evolution can reach regions 
outside the unital region. Here we use the spectator configura- 
tion, S = 0.2, and A = with the GUE model. 
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2 3 4 5 6 7 8 

log2 dimTic 



Figure 4.3: We evaluate the unitality condition for the spectator con- 
figuration {i.e. on one of the qubits). We vary the size of the 
environment and test different times, that would lead approxi- 
mately to different values of concurrence as shown in the inset. 
See text for details. C = 0+ refers to the value a time step im- 
mediately before C = 0. A line with slope —1/2 is included for 
comparison. 
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Figure 4.4: We show D [eq. (4.3 ] which measures a "distance" in the 
CP plane between numerical curves for the RMT model and 
the Wigner curve. In (a) its dependence with the size of the 
environment is studied for two cases. For A = 0.02 (red dots) a 
finite value is approached, whereas for A = 0.14 algebraic decay 
is seen. The black solid line, which is proportional to N^^, is 
meant as a guide to the eye. In (b) we plot D as a function of the 
coupling A, for various values of Ng. In the large environment 
limit, and for A <C 0.1, we have a noticeable deviation from the 
Werner curve. In all cases A = 1. 



environment increases. 

We have established convincingly that the RMT channels are nearly uni- 
tal. For large purities this allows to establish a one to one relation between 
purity and concurrence, as the unital area converges rapidly to the line 
C = P. Another fact remains to be clarified. When/how do the curves 
approach the Werner curve? 

To study this situation in more detail, consider a CP-curve Cnum{P), ob- 
tained from our numerical simulations, and define its "distance" D to the 
Werner curve as 

D= dP|Cnum(P)-Cw(P)|. (4.3) 

It is sensible to compare D with the unital area = 1/18. The behavior 



of D as a function of the size of the system is shown in fig. 4.4 'a). For 
A = 0.14 (black dots), the error goes to zero in an algebraic fashion, at 
least in the range studied. In fact, from a comparison with the black solid 
line we may conclude that the deviation D is inversely proportional to the 
dimension of Tig. By contrast, for A = 0.02 (red dots), D tends to a finite 
value, in line with the assertion that the numerical results converge to a 
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different curve. In fig. 4.4 'b) we plot the error D as a function of A, for 
different dimensions of Hq. The results suggest an exponential decay of 
D with the coupling strength. The simplest dependence of the deviation 
in agreement with these two observations is 

Da=i = 23^/^ + 25iu- (4-4) 



We also plot the curves corresponding to this ansatz in fig. 4.4 Good 
agreement is observed for D <C D^. Notice the exponential decay of D 
with respect to A. One can thus, in an excellent approximation for large 
A, say that for large dimensions the limiting curve is the Werner curve. 
For A = 0, all studied couplings numerical convergence to the Werner 
curve was observed in the large Ne limit. 

In the presence of TRI the CP-curves a similar behavior is observed when- 
ever A > 0. However, in contrast to the GUE case, no saturation of the 
deviation D was observed when A = 0. In the other configurations con- 
sidered (the joint and the separate environment), the behavior is similar. 
In those cases it is the largest (of the two) coupling strength which domi- 
nates the behavior in the CP plane (as well, naturally, as in time). 



4.2 Entanglement decay 

Sufficiently close to P = 1, the above arguments imply a one to one corre- 
spondence between purity and concurrence, which simply reads C ~ P. 
This allows to write an approximate expression for the behavior of con- 
currence as a function of time 

Cir(0 = Plr(0. (4-5) 

using the appropriate linear response result for the purity decay. The 
corresponding expressions for purity decay have been discussed in de- 
tail in the previous chapters. Eq. ( |4.5| l has similar limits of validity as 
the linear response result for the purity. As it follows from a linear re- 
sponse approximation for purity, we call it a linear response expression 
for concurrence decay. 



In those cases where the deviation from the Werner curve [see eq. ( |4.1|] is 
small and where the exponentiated linear response expression eq. (2.21 1 
holds for the average purity, we can write down a phenomenological 
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Figure 4.5: We show numerical simulations for the average concur- 
rence as a function of time in the joint environment configu- 
ration (GUE case). In (a) we consider small couplings Ai = 
A2 = 0.01 which lead to the Gaussian regime for purity decay. 
The red symbols show the result without internal dynamics, 
whereas the black symbols are obtained for fast internal Rabi 
oscillations (Ai = A2 = 10). The theoretical expectation for con- 
currence decay based on eq. ( |4.6| and eq. ( 3.36| l is plotted in 
the corresponding color. In (b) we consider stronger couplings, 
Ai = A2 = 0.1, such that purity decay becomes essentially expo- 
nential (Fermi golden rule regime), while the level splitting have 
been set to Ai = A2 = 0.1 (red symbols). The theoretical expec- 
tation (red solid line) is based on the same expressions as in (a). 
The insets display the corresponding evolution in the CP-plane 
with the same symbols as used in the main graph. In addition, 
the physically allowed region (gray area) and the Werner curve 
(black solid line) are shown. In all cases Ng = 1024 
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formula for concurrence decay, which is valid over the whole range of 
the decay 

Celr(f) = Cw{PELR{t)) ■ (4.6) 



In fig. 4.5 we show random matrix simulations for concurrence decay in 
the joint environment configuration. We consider small couplings Ai = 
A2 = 0.01 which lead to the Gaussian regime for purity decay as well 
as strong couplings Ai = A2 = 0.1 which lead to the Fermi golden rule 
regime. We find good agreement with the prediction of eq. ( |4.6| , except 
for the Gaussian regime when we switch-on a fast internal dynamics in 
both qubits (Ai = A2 = 10) and consequently D is large. See the insets 
in fig. |4.5[ Notice how in both the heuristic formula and the Monte Carlo 
simulations, we observe entanglement sudden death IIYE04I . 

Note that, in the separate environment case, the entanglement of the two 
subsystems defined on the spaces Hi (8) Ti.e and H2 H^' is constant in 
time. For the initial conditions we use, the entanglement stems entirely 
from the two qubits. The concurrence of these two qubits decays because 
the constant entanglement spreads over all constituents of the two large 
subspaces. 

We now investigated pure initial states of the qubits that are not Bell 
states. Consider the GUE case with no internal Hamiltonian, in the spec- 
tator configuration. As discussed in previous sections, thanks to the sym- 
metry of the ensemble, we can write the state as 

It/'ia) = cos0|OO) +sin0|ll). (4.7) 

The initial concurrence is Co = sin 20; for 6 = n/A we obtain initial Bell 
states. Applying the totally depolarizing channel to one or both qubits, 
results in Werner states. For 6 ^ 7r/4 (mod 7z/2), applying the channel 
to a single qubit or to both qubits results in different evolutions in the CP 
plane. Since we are considering the spectator configuration it is sensible 
to apply it to a single qubit. If we parametrize the depolarizing channel 
IINCOOl using the Kraus operators { ^/T^^t, ^7/3 ax, \/ 7/3 dy, v^7/3 a^}, 
we obtain the state 

I ^(3-27)A 1(3-47)7(1 -A)A \ 

1, 



57(1 -A) 

0^0 
V i(3-47)V(l-A)A 1(3-27)(1-A) / 

(4.8) 
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2 4 6 

Figure 4.6: On the top panel the behavior in the CP plane for different 
pure states is shown. We let the initial states have C = 1,2/ 3, or 
1/3. The red curves are the ones described by the states when a 
totally depolarizing channel acts on the first qubit, whereas the 
black dots correspond to Monte Carlo simulations. In the inset 
we see theoretical curves for other initial concurrences. On the 
lower panel in the inset, we see the exponentiated formula for 
purity decay (blue line) and the Monte Carlo averaged purity 
(black dots). On the main panel we see theoretical result (blue 
line) and 30 samples from the numerical ensemble (thin black 
lines). A = 0, A = 0.01, and dimT^e = 256. 
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with Co = 2a/A(1 — A). Concurrence for this state is max{Co(l — 27), 0} 
and purity is 1 + [47^(2 + Co) — 67(2 + Cq)] /9. Inverting gamma in 
terms of purity and the initial concurrence, we can insert it in the for- 
mula for concurrence to obtain the Werner-like relations 




4(1-P) _ 1 ifop > 5_r2 

Cw,Co = Co<|2V 2+^0 2 >5 Co,_ ^^ ^^ 

otherwise 



The comparison between numerics and the ansatz is presented in fig. 4.6 a). 
The agreement is more than satisfactory. We also show in the inset a set 
of curves for more values of Cq. The success of the description allows 
to give a heuristic concurrence decay formula for initial states, not being 
Bell states along the same lines as in the previous paragraphs and using 
eq. ( |4.9| instead of eq. ( |4.1| . A comparison between our heuristic formula 
and Monte Carlo simulation is presented in fig. |4.6| ^b). Note that "sudden 
death" of entanglement IIYE04I persists, though it happens at different 
purities. This is also visible for the individual members of the ensemble. 



,,Co-l , I + 2C0, -1 + Vl-(1 + Cg)(3-6P-Cg 
Cw,Co = S 0' — ^ \ ^ TTT^i 



1+q 

(4.10) 

is the analogous expression for concurrence in terms of purity, when both 
qubits are coupled with similar strength. 



Chapter 5 

An example, the KI chain 



5.1 Introduction 



In this chapter we continue studying decoherence and entanglement de- 
cay using two qubits as the central system. We pursue two main ob- 
jectives. The first one is to test some of the qualitative and quantitative 
conclusions we have derived in the previous chapters, in a deterministic 
model. The second one is to explore dynamical regimes that, up to this 
point, we have ignored: integrable dynamics and intermediate dynamics 
(not integrable but with chaos not fully developed). We shall focus on 
three points. First we wish to study how dynamics of the environment, 
i.e. its integrability or chaoticity influence the behavior of decoherence 
and internal entanglement. Second, we shall analyze to what extent pu- 
rity and concurrence may be related. In particular we want to test the 
reach of the relation discovered in sec. |4] Finally, we wish to relate quan- 
titatively the behavior of purity, in the chaotic regime, with the results 
obtained with RMT in the previous chapters. To implement such a pro- 
gram we need a model with flexible dynamics, which allows efficient 
numerics for large Hilbert spaces. 

The fact that the central system consists of two qubits, makes a spin sys- 
tem an attractive candidate. Indeed the kicked Ising spin chain, a model 
introduced by Prosen [Pro02J, was used to study decay of fidelity and pu- 
rity |PS02I in echo dynamics for integrable, chaotic (more precisely mix- 
ing fProOOl) as well as for intermediate cases. Results for purity in echo 
dynamics can be used for purity decay in standard forward dynamics if 
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the perturbation is chosen as the coupHng between the central system and 



the environment, sec. 2.2 The main advantage of this model is the high 
computational efficiency that allows to perform numerical calculations 
up to twenty qubits IIPSPS06II and more on any good workstation. Yet the 
model in its original form does not allow for variable Ising interactions, 
and is thus not well suited outside of the field of echo dynamics. Exper- 
imental realizations of similar models have been proposed IIBMM05I and 
related studies appeared in IjLSOSp . 

In the present chapter we shall generalize this model allowing arbitrary 
interactions between the spins, and site dependent kicks (this does not 
affect the numerical efficiency of the model, see appendix |B|. We then use 
the generalized model to study the evolution of purity and concurrence of 
central system (consisting of two spins) initially in a Bell state, evolving in 
environments with different dynamical properties. As we want concur- 
rence to be affected solely by the coupling to the environment, we choose 
non-interacting spins for the selected pair. Some preliminary results were 
published in IIPS06L For the environment it is sensible to consider ran- 
dom states to emulate a bath at fairly high temperature. Using unitary 
time evolution of the total system and partial tracing over the environ- 
ment we can then calculate concurrence and purity decay of the selected 
pair, and discuss their behavior. 

To discuss the general model we first recall and discuss a particular case 
of the model which, in fact, will help to understand the properties of the 
environment (sec. |5.2| . Then we state the general model and discuss the 
particular configurations to be studied (sec. |5.3| . Next we proceed to dis- 
cuss the behavior in time of both concurrence and purity for the different 



configurations and different dynamics, sec. 5.4 Afterwards, in sec. 5.5 



we discuss the relation between concurrence and purity. Finally, we com- 



pare purity evolution with the corresponding RMT model in sec. 5.6 



This chapter is largely based on IIPS06I : however, we have learned a lot 
since its publication. Further research and major rewriting were thus re- 
quired to contextualize and discuss the results presented here. We added 
new configurations and discuss in detail the physical differences between 
the different models. The comparison with RMT is now possible as the 
theory is completed. 
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5.2 The kicked Ising spin chain 

We first study the kicked Ising chain (KIC) ||Pro02|: A ring of L spin 1/2 
particles which interact with their nearest neighbors via a homogeneous 
Ising interaction of dimensionless strength /, and are periodically kicked 
with a dimensionless homogeneous magnetic field h.The Hamiltonian is 
thus 

W = /E^/t^^l+^l(OE^-^; (5.1) 

where 5i{t) = Ylnez^i^ ~ represents an infinite train of Dirac delta 
functions with primitive period one; a^'^'^ are the Pauli matrices of par- 
ticle i and aj = {a-^ , crj , aj) . We must also impose periodic conditions in 
order to close the ring: (7^+1 = ^^i- During the free evolution, i.e. between 
the kicks, the system evolves with the unitary propagator 

L7ising(/) = exp (^-ij E cr/cT/^i^ , (5.2) 
and the action of the kick is described by the unitary operator 

Ukick(^) = exp -i E ^ • '^i^ ■ (5.3) 
The Floquet operator for one period is thus 

UkI = LJlsing(/)lJkick(^)- (5.4) 



A detailed discussion of the symmetries of the system can be found in 
||PP07|. We summarize the main observations here. Due to the homo- 
geneity of the kick b and the nearest neighbor interaction /, the system 
has a rotational symmetry iij — > o'j+i- The dimension of each of the invari- 
ant subspaces (Hk with A: = 1, . . . , L) is close to 2^/L for large L. Each of 
the invariant subspaces has an anti-unitary symmetry which, if we choose 
a base in which b ■ (ij is real, is complex conjugation (i.e. not conventional 
time reversal symmetry). Since each minimal invariant subspace has an 
internal anti-unitary symmetry, it is reasonable to compare with the GOE 
or the COE to understand the properties of the system. 
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Figure 5.1: The nearest-neighbor spacing of the homogeneous spin 
ring (as a black curve) with a magnetic kick b = (0.8,1.4,0). 
We compare with the Brody distribution P{co,s). Its limiting 
values CO = 1 and a; = yield GOE and Poissonian statistics re- 
spectively. Intermediate values suggest intermediate dynamics. 
Here we use 15 qubits. 

As mentioned above the system is very rich regarding its djmamics. Var- 
ious parameter combinations yield integrable (meaning analytically solv- 
able) dynamics. In particular setting the magnetic kick parallel or per- 
pendicular to the preferred direction of the Ising interaction implies inte- 
grability. However, using b parallel to the Ising interaction yields trivial 
integrability (all aj commute with the Hamiltonian), so we shall use b 

perpendicular to the z axis. In particular we shall set b = (0, 1.53,0) since 
for this parameter some oscillations that we shall observe have a bigger 
amplitude and period. In order to obtain chaotic dynamics one cannot 
just choose parameters outside the set that yield integrable dynamics. 
One must be "far enough" from this set. We verified that for / = 1, 
b = (1.4,1.4,0) (used in this chapter) and / = 0.7, b = (0.9,0.9,0) (used 
in the following one) the system is chaotic (see |PSPS06l IPP07I ): many 
statistical tests like the nearest-neighbor spacing, the form factor, and the 
skewness, were applied. Except at extremely large energies, no signifi- 
cant deviations from RMT were observed. In order to obtain intermediate 
dynamics we wanted to be far enough from both integrability and com- 
pletely mixing dynamics to have an appreciable compromise regarding 
spectral statistics. For / = 1 and b = (0.8,1.4,0) we indeed find interme- 
diate spectral statistics. Their nearest-neighbor spacing is characterized 
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by a Brody parameter pBro73] of 0.33, see fig. 5.1 The nearest-neighbor 
spacing distribution in the Floquet spectrum shows a marked level repul- 
sion, but it is still quite far from the distribution we expect and get for 
the chaotic case. 



5.3 A generalized kicked Ising model 

The Hamiltonian of the generalized kicked Ising model is 

E /M^/^ + '^i(OEV^r (5-5) 

j>k^i j=i 

The model thus consists of a set of L spin 1/2 particles coupled to all 
other spins by an Ising interaction (first term) and periodically kicked by 
a site dependent tilted magnetic field (second term). 

Our model differs from the one stated in eq. by the fact that the cou- 
pling Jjjt^ is between any pair of particles and has an arbitrary strength. 
This freedom allows to build a variety of models with different topolo- 
gies. Each of them will correspond to a well defined physical picture, and 
thus we show how the flexibility of the model can be exploited to study 
diverse physical situations in a common framework. We also allow site 
independent kicks to be able to have (or not) internal dynamics in our 
central system. 

The central system is composed of two spins, say spins "1" and "2". 
We are left with = L — 2 spins which are going to be considered as 
the environment. We furthermore diminish the number of parameters 
by setting &i = &2 = arid bs = ■ ■ ■ = bi = bg. The fact that we 
keep the kick in the central system can represent local operations made 
by the "owners" of each of the qubits, and will not affect the values of 
concurrence and purity. We shall furthermore require weak coupling of 
the central system to the environment. Thus 7,4 and }i^2 must be much 
smaller than the typical Ising interaction within the environment. We 
set 72,1 = in order to prevent any interaction between the spins in the 
central system. 



We propose in particular six models. They are schematically drawn in 
fig. 5.2 The open circles represent qubits 1 and 2 (the central system) 
whereas the filled circles represent the other qubits (environment). The 
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(a) (b)g (c)^ 




Figure 5.2: Different configurations of the coupling of the central sys- 
tem (represented as open circles) to the spins that act as envi- 
ronment (represented as filled circles). Thick/thin lines repre- 
sent /ce//e respectively. The description for each model and its 
features is given in the main text. 



thin lines represent a weak Ising interaction /^g of a particle in the central 
system with one in the environment. The thick lines represent an interac- 
tion /e among particles in the bath; in all numerical results presented we 
fix /e = 1. Finally we note that the dimension of the Hilbert space of the 
environment is Ng = 2'?«. Writing the Hamiltonian for each of the models 
we propose is now a trivial task. As an example we write the Hamilto- 



nian of the model (a) in fig. 5.2 Rewriting the Hamiltonian in terms of 
central system, environment and interaction as H = Hc + Hg + Hce the 
parts are given by 

2 

Hc = 3i{t)J2bc-^j, (5.6) 

He = Je E ""M+l + '^l (f ) E • (5.7) 

j=3 j=3 

Hce = ]>ICTI. (5.8) 



Each of the models has a characteristic Heisenberg time Th. For a sys- 
tem with no symmetries the Heisenberg time is defined simply as In/ d, 
where d is the mean level spacing. If the system has symmetries, for 
each sector we have a Heisenberg time defined again as 2n/ d, but d is 
the mean level spacing in each sector. We shall stress the Heisenberg 
time of each model, as our RMT models show a strong dependence on 
it [e. g. eq. ( 3.23| ]. We also use a normalized coupling /ce, which gives 
a model independent coupling strength. Finally, we want to mention the 
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corresponding configuration (spectator, joint environment or separate en- 



vironment), as in sec. 3.1 



Please refer to fig. 5.2 while reading the descriptions. 



• Model (a) is in the spectator configuration. A single qubit is cou- 
pled to the extreme of an open chain. Due to the breaking of the 
symmetry of the coupling, the Heisenberg time of the environment 
is th ~ Ne. We use the normalized coupling Jce = Jce- 

• The second model is a variation of the previous one. In this the 
second qubit is coupled to the other extreme of the chain. It is an 
example of the joint environment configuration. Again Th ~ Ng, 
however to have the same effective coupling we must scale Jce = 

Vlj'ce- 

• Model (c) is the last of this family of models in which the environ- 
ment is an open chain. In this case Th ~ \/Ne ^ Ne. As we have 
two coupling spots, we must scale Jce = V^Jc^. This is a good exam- 
ple of a separate environment Hamiltonian. This, and the previous 
two models, were analyzed in IPS06 I . 

• Model (d) is again in the spectator configuration, but here we have 
a ring (thus the translational symmetry is not broken) and the cou- 
pling is symmetric. The Heisenberg time is then related to each 
symmetry sector separately and is approximately equal to Ng/cje- 
The coupling must be normalized to Jce = ^/^eJ'ce- 

• The next model differs from the previous one only inasmuch that 
the coupling is not symmetric. This is enough to mix the invariant 
subspaces and thus the Heisenberg time is ~ Ng. There is a single 
coupling and thus it is enough to set Jce = Jce- 

• Model (f) is another example of a separate environment config- 
uration. It has the smallest Heisenberg time, since the coupling 
does not break the symmetry of the ring. Thus, it is approximately 
V^Ne/ qe- The coupling is normalized to Jce = ^/'Je /ce- 
lt is pertinent to mention, that the degree of chaoticity was obtained for 
cyclic chains, while the chains representing some of our environments are 
open. Yet, for large number of spins, this is irrelevant. 
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The non-unitary evolution of the central system alone is calculated along 
the same lines as for the RMT models. We perform a unitary evolution 
of the whole system (yielding state \ip{t))) and then we partial trace over 
the environment. That is 

Pit) =tr, \tp{t)){ip{t)\. (5.9) 
We shall consider an initial condition 

\xp{t = 0)) = \ipBea)(^\lpenv), (5.10) 

i.e. a general Bell state [any state such that C(|i/?Beii) ('/'Belli) = 1] riot 
entangled with the environment, which is in a pure random state [for 
models (a), (b), (d), and (e)] or in a tensor product of pure random states 
in each section of the environment [for models (c) and (f)]. 



5.4 The evolution of concurrence and purity 

We now present the results of our numerical calculations of both concur- 
rence and purity of the selected pair of spins as a function of time. We 
first concentrate on the short time behavior. 



Fig. 5.3 shows the time evolution of purity of the selected pair of spins for 
the initial state consisting of a Bell pair in the central system and the ap- 
propriate random state in the environment. We choose the environment 
to be chaotic in one case and integrable in the other, using the correspond- 



ing parameters mentioned in sec. 5.2 The chaotic environment leads to 



linear decay for some time interval, while the integrable environment 
leads to a quadratic decay. The interval in which the linear decay can be 
observed depends on two things: the configuration and the dimension 
of the environment. Indeed, for configuration (f) we have the shortest 
linear decay but we also have the shortest Heisenberg time. The reader 



is encouraged to check table 5.4 to understand better our numerical re- 
sults. Recall that in all our derived formulae, we predicted linear decay 
before Heisenberg time. The situation for the integrable case is different. 
We observed a marked quadratic decay, independent of the configuration 
and the size of the environment. We believe that one should be able to 
explain the results (for the integrable situation) with exact analytic cal- 
culations, however the aim of this chapter is not to explain the details, 
but to illustrate the qualitative differences between having a chaotic and 
regular environment. 



5.4. The evolution of concurrence and purity 
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Figure 5.3: These figures show the evolution of purity for a sir\gle 
initial condition per configuration, for chaotic (points) and inte- 
grable (thin lines) dynamics. All different configurations from 



fig. 5.2 are tested (we keep the corresponding color coding). 
In the chaotic case, purity shows a linear dependence on time 
whereas in the integrable one, it exhibits a quadratic decay. We 
also plot the line 1 — 3}^t and the parabola 1 — with thick 
black curves, as guides. We vary the size of the environment 
to understand the scaling with q^. The corresponding plots for 
concurrence are indistinguishable from the ones presented here. 
We have set k = 0.005 and be = h- 
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(a) 


(b) 


(c) 


(d)B 


(e)B 


(f)" 


qe = 12 


4096 


4096 


64 


341.3 


4096 


10.7 


qe = 14 


16348 


16348 


128 


1170.3 


16348 


18.3 


qe = 16 


65536 


65536 


256 


4096 


65536 


32 


qe = 18 


262144 


262144 


512 


14563 


262144 


56.9 



Table 5.1: Numerical values of the Heisenberg times for the different 
configurations of the environment, and different sizes of it. 



We now proceed to look at the long-time evolution of concurrence and 



purity again for all configurations. In Fig. 5.4 we show both concurrence 



and purity for integrable, intermediate and chaotic situations. 

For the integrable case we see that for some configurations [(a), (b), and 
(c)] after the initial quadratic decay there is a full revival of both concur- 
rence and purity. However the autocorrelation function of the full system 
(not shown) drops fast and reaches negligible values even after the first 
few kicks. It is also worthwhile to note that the oscillations of purity 
imply increase in some cross-correlation function fPSZ03|. For the other 
configurations there is no such revival, and the decay seems Gaussian; 
this is in agreement with the quadratic decay observed for large purities. 
Notice that in all the configurations in which there is a revival the cou- 
pling is to one end of a chain. Notice how concurrence drops to exactly 
after a finite time. 

For the chaotic case we observe a monotonic decay of both purity and 
concurrence. The speed at which they decay, though initially the same, 
is quite different for the different configurations. This can be explained 
with the huge differences in the Heisenberg times of the environment. 

For the mixed case, we observe a compromise among the behavior in the 
integrable case and the chaotic one. Interestingly in all cases we observed 
that concurrence drops identically to zero after some time, in contrast 
with purity which does no reach its minimum. 



5.5. A relation between concurrence and purity 
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Figure 5.4: Long-time behavior of purity (left) and concurrence 
(right) for an arbitrary initial condition in integrable (top), non- 
ergodic and non-integrable (middle), and fully chaotic (bottom) 
regimes. For the integrable regimes we observe recurrences for 
the open environment configurations and fast, possibly Gaus- 
sian, decay in the others. For the intermediate case we observe 
in some instances the superposition of oscillations together with 
the exponential decay. For the chaotic case we observe expo- 
nential decay together with a marked dependence on the kind 
of environment. We have set = 16, }c = 0.02, and be = bg. 



5.5 A relation between concurrence and purity 



The rather similar pictures emerging for concurrence and purity suggest 
that a simple relation between the two might emerge. Indeed, from the 
results of the previous chapter this could be expected for the chaotic case. 
Unfortunately we are not able to access such a big range of values of the 
coupling as in the previous section; too strong coupling leads to total con- 
currence decay after a single period (remember that we are dealing with 
Floquet operators). Too weak couplings imply too long runs as we have to 
integrate one step at a time. On the other hand, we are going to be able to 
test the relation not only on chaotic situations but also in other dynamical 
regimes such as integrable environments and intermediate environments. 
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Figure 5.5: The evolution of a random initial condition, in the (C, P) 
plane, for different dynamical regimes: integrable, intermedi- 
ate, and chaotic, for two different coupling strengths. All con- 
figurations are tested, the color indicates the configuration, as 
in fig. 5.3 In the thermodynamical limit all curves seem to fall 
in the unital region. For smaller couplings there is a tendency 
to stick to the Werner curve. For stronger couplings, in the in- 
tegrable and intermediate regime, they seem to fill the unital 
area. We set here = 16, and be = bg. The gray region and the 
black thick line have the same meaning as in fig. 4.1 



5.5. A relation between concurrence and purity 



65 



We plot in fig. 5.5 concurrence against purity in all the dynamical regimes 
studied so far. We do so for two strengths of the perturbation, namely 
}c = 0.1 and 0.01. A wide variety of behaviors are observed. One thing 
that stands out is that all curves lie in the unital region. Indeed for / = 
0.1, in the integrable and intermediate regime we observe a filling of the 
unital region [recall it from fig. |4.1| . This is not a finite dimension effect: 
we verified that the behavior remains for increasing dimension of the 
environment in both cases. For small dimensions, and in configuration 
(b) one observes that the curves do not lie in the unital region and explore 
a larger area. This can be explained: one has a common environment 
which, in fact, has memory (due to its small dimension). An example of 
this behavior can be found in figure 5 of l|PS06| . 

For the other cases, where there is no filling of the unital region, one can 
suspect self averaging. This can be tested calculating the standard de- 
viation of purity dp , one time step before concurrence is zero. This is 



equivalent to examining the width of the coloured curves in fig. 4.1 We 
numerically observed that for some set of parameters, ap ck l/Ng whereas 
for other it approaches a finite value. Though a tendency to accumulate 
towards the Werner curve is clear from the pictures presented, in some 
cases it is not more than a tendency. Indeed, even in the chaotic case and 
when having self averaging, the limiting curve not always crosses the line 
C = 0atP = l/3 (like the Werner curve), but a slighter bigger value. 
Again our numerics show that this is not a finite size effect. However, all 
the limiting curves are similar to the Werner curve, and that still can be 
used for approximations, as long as we recall it is only an approximation. 
Though all six configurations are physically quite different and the indi- 
vidual behavior of purity and concurrence is quite affected, the relation 
between the two is entirely robust. 

Cases (c) and (f) lead to an instructive lesson: Here we have two uncou- 
pled environments, and we start with a pure state in each of these. The 
purity of the uncoupled subsystems will remain unchanged, but the pu- 
rity and concurrence of the initial Bell pair will decay. Thus one might 
consider seeing a paradox, but this is not the case; the entanglement of 
the pair is simply spread over all the system with time. 

Up to this point our results are in full qualitative agreement with the ones 
found in previous chapters. This leads us to make a tighter comparison 
with the RMT models developed before. 
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5.6 Comparison with RMT predicted behavior 



Our KI model has some parameter ranges in which we can consider the 
system to be "quantum chaotic". A quantitative comparison of purity 
decay for the dynamical system and the RMT model is sound. Statistical 
tests and analytical arguments, indicate that one can consider the homo- 
geneous KI as a typical member of the GOE. Thus we shall compare with 
the results obtained in sec. |2.5l and sec. 13.2.21 We first list some differences 
between the two models. 

• The most serious difference is the structure of the coupling. Relations 
of the form eqs. ( |E.4| and ( |E.5| used to define the coupling in the RMT 
model, generically yield non-separable operators, whereas the coupling 
for the KI model is indeed separable and of the form (7^ (g) Ve- Dynamics 
involving separable operators have special properties. We must not be 
surprised that this particular structure of the coupling affects the behav- 
ior. 

• The theory developed in earlier chapters is for an ensemble of systems. 
Here we are dealing with a single system. Self averaging takes care of 
this issue. Numerical tests indicate that as the size of the environment 
increases, the behavior for a typical Hamiltonian approaches rapidly the 
behavior of the ensemble average (see fig. |2.3| |. We can safely ignore 
this problem as long as we go to high enough dimensions. However, 
for very detailed discussions this must be taken into consideration (as in 
IIPSPS06II ). 

• We are using environments with either exact rotational symmetries [(d), 
(e), and (f)] or with some approximate translational symmetry [(a), (b), 
and (c)]. But we have not developed a theory for chaotic environments 
with discrete unitary symmetries. 

One might be tempted to assume that the overall result is an average of 
the results in the different sectors, as happens with fidelity /: Given an 
echo operator M with a discrete unitary symmetry, and a set of states 
{|j/'s)} living in the different symmetry sectors. 



S s' s 
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For purity the situation is different, 

V s,s' / s 

Let ^1 = ^2 = 1/ V2, \'^s=i) = |00), and \ips=2) = In this case purity 
of the reduced system is minimal (1/2) whereas the purity of the reduced 
state of each of the components is maximal (1). 

However, a quantity that is affected immediately when superposing spec- 
tra is the form factor. The form factor of a superposition of many, say 
GOE, spectra approaches that of a Poissonian one. Thus, a simplified 
way to incorporate this effect is to replace in our formulae the GOE form 
factor &2^^ {t) by the one of a random spectra &2(0 = 0. 

• During the first chapters we assumed time independent Hamiltonians. 
The KI Hamiltonian ( |5.5| l depends explicitly on time. However, due to 
its periodic nature, one can construct its corresponding Floquet opera- 
tor. The circular ensembles |Dys62| are used to model such operators 



IIHaaO ll. Circular and Gaussian ensembles share all the properties re- 
garding its eigenvectors. The eigenphases of the former and the eigenval- 
ues of the later share most properties but a precise mathematical relation 
among them is cumbersome to state. We can naively consider the circular 
ensembles as an exponentiated (and approximated) version of the Gaus- 
sian ensembles. We thus do not need to modify our theory to take into 
account this issue. 

• With eq. ( |E.4| we are assuming certain normalization of the coupling. 
One can adjust the coupling constant A so any V has the given normal- 



ization [see eq. (2.5 1]. In our case, we can expect that for models (a) and 
(e) the decay, at least for short times, does not depend on the number 
of qubits. This can be achieved substituting A ocjce/ \/^ with a. an 
Ne independent number. The leading term (with respect to time) is now 
independent of Th and Ng. 

For simplicity we eliminate the internal Hamiltonian in the dynamical 
model (He = 0) to use the RMT formulae with A = 0. We use config- 
uration (d) for the comparison to eliminate possible border effects. Its 
Heisenberg time is moderately small and thus we can access easily time 
regimes comparable with the Heisenberg time. The effective size of the 
environment is still large (2'''^) [compared with separate environment con- 
figurations (v^2^)] so the large dimension limit is achieved faster. 
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Figure 5.6: We plot, with blue thin lines, the evolution of purity for 
10 initial condition in model (d). We vary the number of qubits 
vertically and the time ranges horizontally. The formula de- 
rived from our RMT model eq. ( |5.12 1 is plotted as a thin black 



line. The same formula, but assuming no correlations within 
the spectra, is plotted as a dashed line. We use /^g = 0.0005 
to be in the linear response regime for all plots. See text for a 
complete description. 

Our adapted RMT formula (using Bell states as initial conditions) reads 



The numerical results are summarized in fig. 5.6 There we find 9 plots 



for the evolution of purity in the KI chain. Each blue line corresponds to a 
different realization of the initial condition eq. ( 5.10| . We include formula 
(5.12 1 with and without the correlation term ^B^^t) as solid and dashed 
black lines respectively. We increase the size of the environment (from 
bottom to top)to appreciate how the large dimension limit is approached. 



5.6. Comparison with RMT predicted behavior 



69 



We also vary the plot range to examine different time regimes of the 
problem. In the plots, two time scales are available. The lower scale is the 
number of time steps (number of applications of the Floquet operator) 
while the upper one is time in units of the Heisenberg time. 

We first encounter the Zeno regime, in which the discreteness of the spec- 
trum induces deviations of the RMT formulae; for a finite spectrum, what 
we take as a 3 function in eq. ( 2.13| , has a finite width and height, caus- 



ing an initial quadratic decay. Indeed we can see on the left most plots 
that the deviation lasts a few kicks. The effect is more visible for = 12 
and 14. The next regime is the linear one, sometimes called Fermi golden 
rule regime. We observe linear decay as predicted by RMT whenever 
1 <C i <C Th. We use this regime to obtain the fitting parameter a = 0.21. 
Next we observe a crossover in which deviations from the linear behav- 
ior are evident (starting at t ^ Th/4). For = 10 this regime overlaps 
slightly with the Zeno regime. After the crossover we find the Gaussian 
regime, in which a quadratic decay of purity is manifest. The clearest 
picture of this behavior is seen for q^ = 10 where the Heisenberg time is 
the smallest. 

As we increase the size of the environment we observe better agreement 
between the dynamical model and the RMT formula with no correlation. 
For q^ = 14 the agreement is virtually perfect for all time regimes beyond 
the Zeno regime. Comparing the behavior for different environments, we 
can observe the dramatic effect of the Heisenberg time on decoherence: 
The denser the spectrum, the later we arrive to the Gaussian regime. En- 
vironments with tighter spectra cause less damage to the central system. 
Other configurations also show qualitative agreement with the RMT for- 
mula developed before. 



The results of sec. 2.5 predict that for arbitrary separable conditions, on 
a TRI system (or with any other anti-unitary symmetry), we have a fi- 
nite spreading of purity, even in the large dimension limit. This effect 
becomes important after the Heisenberg time. We argued why we expect 
self averaging (with respect to the initial conditions) for Bell states and 
indeed observed it in the large dimension limit. 

We tested these two findings in our KI models. As expected, we found no 
self averaging for separable ones. However, we observed the same effect 
before and after the Heisenberg time. This anomaly can be due to the 
separable structure of the coupling, which implies a preferred base within 
the central system. From that point of view, this irregularity is analogous 
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Figure 5.7: Evolution of 25 random initial conditions for configura- 
tion (d). In the top panel we prepare separable initial condi- 
tions, whereas in the lower panel entangled ones. We observe 
no self-averaging both for long times and short times (in the 
insets). This might be due to the separable structure of the cou- 
pling, /ce = 0.0005 and qe = 12. 

to the effect of an internal Hamiltonian. Interestingly, we again observe 
that for entangled states we do have self averaging. This suggests that the 
"transport" of symmetries via entanglement is a general property that is 
not limited to the topics discussed during this thesis. 



Chapter 6 

Quantum memories 



Divide et impera. 



Physical devices capable of storing faithfully a quantum state, i.e. quan- 
tum memories, are crucial for any quantum information task. While dif- 
ferent types of systems have been considered, most of the effort is concen- 
trated in manipulating qubits as they are essential for most quantum in- 
formation tasks I.NCOO I. For a general quantum memory (QM) it is neces- 
sary to record, store, and retrieve an arbitrary state. For quantum commu- 
nication and quantum computation initialization in the base state, single 
qubit manipulations, a two qubit gate (e.g. control-not), and single qubit 
measurements are sufficient. For arbitrary quantum memories recording 
and retrieving the state is still an experimental challenge but initializing 
a qubit register and individual measurement is mostly mastered. Faithful 
realization of two qubit gates together with better isolation from the envi- 
rorunent are the remaining obstacles for achieving fully operational quan- 
tum technology. However the huge effort done in the community has 
born some fruit lJSC+041 ICdRF+05l iCMJ+OSl lEAM+051 lPF02l IPFL+Olt. 



Understanding decoherence of one, two, and n qubits has been of interest 
for some time. In IIFFP04I the authors calculate decoherence of a QM 
for a spin-boson model using a measure designed for their purposes. 
In the limit of slight decoherence they obtain various results including 
additivity. Experimental results are also available IIKS041 IKS06I . In the 
previous chapters we have studied one and two qubit decoherence using 
both a random matrix formulation and a specific dynamical model. One 
would wish to encompass some of this progress in a general picture. 
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In the present chapter we address this problem using a standard mea- 
sure of decoherence, namely purity. We obtain analytic expressions for 
the decoherence of a QM during the storage time. Specifically we dis- 
cuss a QM composed of a set of individual qubits interacting with some 
environment. The expressions given are based on previous knowledge 
of the decoherence of a single qubit entangled with some non-interacting 
spectator. Their validity is limited to small decoherence, i.e. large pu- 
rity of the QM. Note that the latter is not a significant restriction, within 
the context of quantum information processing, due to the high fidelity 
requirements of quantum error correction codes. We again assume that 
the entire system is subject to unitary time-evolution, and that decoher- 
ence comes about by entanglement between the central system (CS) and 
the environment. Spurious interactions inside the central system are ne- 
glected. 

A further and critical assumption is the independence of the coupling 
of different qubits with the environment in the interaction picture. This is 
justified if the couplings are already independent in the Schrodinger pic- 
ture or if we have rapidly decaying correlations due to mixing properties 
of the environment |PS02|. Physically, the first would be more likely if 
we talk about qubits realized in different systems or degrees of freedom, 
while the second seems plausible for many typical environments. 

The central result is a decomposition of the decoherence of the full QM, 
coupled to a single or several environments into a sum of terms. Each 
of these describes the decoherence of a single qubit in a "spectator con- 
figuration" eq. ( |1.2| . Recall, the CS consists of two non-interacting parts, 
one (the qubit) interacting with the environment and the other (the rest 
of the QM) not. This configuration is non-trivial if the two parts of the 
CS are entangled. Apart from the above assumptions, this result does not 
depend on any particular property of the environment or coupling. Thus, 
it can be applied to a variety of models. 

The general relation is obtained in linear response approximation and 
leads to explicit analytic results if the spectral correlations of the envi- 



ronment are known, see sec. 6.1 We test successfully the results in our 
random matrix model (sec. 6.2 1. Finally, we perform numerical simula- 
tions for four qubits, interacting with the kicked Ising spin chain as an 
environment, in sec. |6.3| 



6.1. The calculation 
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6.1 The calculation 



We now outline the calculation for n qubits. Since it just represents a 
variation in the indices of some calculations presented in appendix |A| we 
skip some details and go through it pretty fast. 

The central system (our QM) is composed of n qubits. Thus, its Hilbert 
space is T^qm = (8)"=i 'Hj, where Hi are the Hilbert spaces of the qubits. 
The Hilbert space of the environment is again denoted by T^e- The Hamil- 
tonian reads 

H = Hq^ + He + kV, Ay = A,y('). (6.1) 

Here, Hqm = Y^^i H„ where H, acts on Hi, whereas He describes the 
dynamics of the environment, and XV the coupling of the qubits to the 
environment. The strength of the coupling of qubit i is controlled by 
the parameter A,, while A = max{|A;|}. The (possibly time-dependent) 
Hamiltonian gives rise to the unitary evolution operator U\{t). 

We consider two different settings. In the first one V''' acts on the space 
Ti-i (8) Ti-Q, i-e. all qubits interact with a single environment called joint 
environment. In the second one each qubit interacts with a separate envi- 
ronment. Th us t he environment is split into n parts. He = 0[Li T^e,i ^rid 



in eq. ( |6.1| , acts only on Hi (8) He,i- The first case would be typical 
for a quantum computer, where all qubits are close to each other, while 
the second would apply to a non-local quantum network. Both configu- 
rations are defined in perfect analogy with the ones, with corresponding 
names, considered in sec. |3.1| 



Again we choose the initial state to be the product of a pure state of the 
central system (the QM) and a pure state of the environment 

\^0) = \^qm)\tpe). \lpqm) G Wqm, \rpe) S He- (6.2) 

In the separate environment configuration we furthermore assume that 
\^e) = ^i \^e,i) with \tpe,i) ^ 'Hgj, which corresponds to the absence of 
quantum correlations among the different environments. 

Purity evolves as P{t) = trpl^{t), with pqm{t) = tr^ Ux{t)\tpo) {xpo\Ul{t). 
To calculate purity (or any other quantity that depends solely on the 
Schmidt coefficients) we replace the forward time evolution Ux{t) by the 
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echo operator M(f) = Uo{t) U\{—t) where Uo{t) gives the evolution with- 
out coupling, recall sec. |2.2[ Using Born expansion for the echo operator 
and defining \Vt = Uo{t)A-VUQ{t) (the coupling at time t in the interac- 
tion picture) purity is then given by 

P{t) = 1 - 2A2 dr dr'Re A{t, t') + O (^A^) , (6.3) 

with 



A{t,t') = p[VrVr'Qo (8) ^o] - pIVt'QqVt (S> Qo] 

+ p[VrQO O Vr'Qo] " P[Vt'QO ® QoVr], (6.4) 

and Qo = \xpo){tpo\ [recall that p[pi pi] = tr(tre|Oi trepi)]- Note that 
the linear terms vanish identically after tracing out the environment. 
Considering the form of the coupling in eq. ( |6.1| l we can decompose 

A^A{t,t') = A,AyA(''^)(T,T') with 

a('V)(t,t') = p[y«y^/'(?o ^ Qo] - p[vi:\ovP Qi] 

+ p[#\o ^ V^pQo] - p[V^^Qo ^ QoVP]. (6.5) 



y/') = LJ+(i)y(OLJo(f) in analogy with Vf Equation (|63 IS given as a 



double sum in the indices of the qubits. In a diagonal approximation 

(i = i), P{t) is expressed in terms of the purities Pgp [t) which correspond 
to the purity decay of the CS in a spectator configuration where only 
qubit i is interacting with the environment. Purity then reads as 

P{t) = (6.6) 

i=i 

P§ {t) = 1 - 2Af 1^ dr I* dr'M'''^ {r, r') + O (a^ 

This is our central result of this chapter. The diagonal approximation 
is justified in two situations: First if the couplings V^'^ of the individ- 
ual qubits are independent from the outset, as would be typical for the 
separate environment configuration or for the random matrix model of 
decoherence. Second, if the couplings in the interaction picture become 
independent due to mixing properties of the environment, as would be 
typical for a "quantum chaotic" environment. 



6.2. A RMT example 
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6.2 A RMT example 

To illustrate the case of independent couplings mentioned above, random 
matrix theory provides a handy example. Such models were discussed in 
previous chapters and in llPS07llPGS07llGS02llAV07l . 

In chapter|3]purity decay was computed in linear response approximation 
for two qubits, one of them being the spectator. For the sake of simplicity 
we choose the joint environment configuration, no internal dynamics for 
the qubits, and He and V^^'^ as typical members of the Gaussian unitary 
ensemble. Purity is then given by 

Pffit) = l-A?(2-pi)/(0; (6.7) 
/(f) = f max{f,TH} + ^ min{t,TH}l (6.8) 

Here, Th is the Heisenberg time of the environment and pi is the initial 
purity of the first qubit alone, which measures its entanglement with the 
rest of the QM, see eq. ( 3.23| . As we required a priori independence of 



the couplings, we can now insert eq. ( |6.7| in eq. ( |6.6| to obtain the simple 
expression 

P{t) = l-f{t)t^ii2-P^ (6.9) 

/=i 

where p, is the initial purity of qubit z. In the presence of internal dy- 
namics the spectator result is also known and can be inserted. As an 
example, we apply the above equation to an initial GHZ state (( |0 • • • 0) + 
\1---1))/V2). Then all Pi = 1/2 and we obtain P(i) = 1 - (3/2)/(i) A^. 

For a W state (|10 . . . 0) + |01 . . . 0) H + |00 . . . purity for each 

qubit is p,- = {n^ — 2n + 2)/n^ 1, in the large n limit, and purity decays 
asP{t) = l-f{t)Ei^l 



6.3 A dynamical model 



We now use the homogeneous kicked spin (sec. \5.2\ as an environment 
for n qubits. We consider its chaotic and integrable regimes. We shall 
thus use a particular realization of eq. (5.5 1 in the spirit of chapter |5] 



The main assumption in eq. ( |6.6| is the fast decay of correlations for 
couplings of different qubits to the environment. For the random ma- 
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Figure 6.1: Purity decay for a GHZ state is shown for the environ- 
ment described by eq. ( 6.10| |. The couplings are: all qubits to 
one spin (red triangles), nearby spins (blue stars), and maxi- 
mally separated spins (green squares). The theoretical result 
(thick line) is calculated in terms of the Pj^^^-j (0 (thin lines). The 
main figure corresponds to the case of a mixing environment 
and the inset for an integrable one. 
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trix model discussed above this is trivially fulfilled. Yet, integrable en- 
vironments are commonly used [CL83| and one may wonder whether 
eq. ( |6.6| can hold in such a context. We shall therefore study a dynamical 
model where a few qubits are coupled to an environment represented 
by a kicked Ising spin chain using identical coupling operators for all 
qubits. In this model, the variation of the angle of the external kick- 
ing field allows the transition from a "quantum chaotic" to an integrable 
Hamiltonian for the environment |Pro02[ IPP07| (a). For convenience we 
separate it in the environment and the central system. The Hamiltonian 
of the environment, namely the homogeneous kicked Ising spin chain, 
recall, is given by 



H(-=) = J2 c^^'V^'+'^ + Sit) J2 b ■ <7^''^ (6.10) 

/ = ! = 



where S{t) = J^nez ^(^ ~ i^-^- time is measured in units of the kick pe- 
riod), L the number of spins in the environment, cr^*^'') = {ax'^''\ cry^''\ (Tz^''^ ) 
the Pauli matrices of spin z, and b the dimensionless magnetic field with 
which the chain is kicked {h = 1). We close the ring requiring p-^-^'^^ = 
cr{e,o) Tj^g Hamiltonian of the QM is 

^(qm) ^ S{t)Y2b- (7(1"^''), (6.11) 
i 

where a^V^''^ is defined similarly as for the environment. The coupling is 
given by 

AV = A^^^'^'V^^''^ (6.12) 

The 7,'s define the positions where the qubits of the QM are coupled to 
the spin chain. Equivalently we could use kicked Ising couplings and a 
time-independent field MPSWOll . 

To implement a chaotic environment, we use here b = (0.9,0.9,0). We 
use a ring consisting of L = 12 spins for the environment and 4 addi- 
tional spins for the qubits of the QM. The coupling strength is A = 0.005. 



In Fig. 6.1 we study purity decay when all four qubits are coupled to 
the same spin, to neighboring spins, and to maximally separated spins, 
; = (0,3,6,9). The initial state is the product of a GHZ state in the QM 
and a random pure state in the environment. We compare the results 



with eq. (6.6 1 (thick line) obtained from simulations of the spectator con- 



figuration (thin solid line). Coupling the spectator to different positions 
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in the chain yields near identical results for P{t), so we can see only one 
line. The figure demonstrates the validity of eq. ( |6.6| for well separated 
and hence independent couplings. For coupling to neighboring spins de- 
cay is slightly faster while the sum rule does not hold if all qubits are 
coupled to the same spin. 

A similar calculation for integrable environments with b = (0,1.53,0) 
yields faster purity decay than in the mixing case as expected from gen- 
eral considerations in [GPSZ06J. Nevertheless the sum rule is again well 
fulfilled except if we couple all qubits to the same spin. This leads us to 
check the behavior of the correlation function 

P[V^''^VIPQ0 ^ Qo] = i^olV^'^V^Pl^o). (6.13) 



6.2 



This is the simplest and usually largest term in A^''i\ eq. (6.5 1. Figure 
shows this quantity for chaotic (i.e. mixing) and integrable environments 
in the first and second row, respectively. The first column shows the 
autocorrelation function (/ = j). The second and third columns give the 
cross correlation function (z ^ j) when the qubits are coupled to the 
same or opposite spins, respectively. In the latter case, correlations are 
always small thus showing that also in integrable situations our condition 



can be met. Non-vanishing correlations, as in fig. |6.2|^b) and fig. 6.2 'e) 



lead to deviations from the sum rule. The pronounced structure for the 



integrable environment, fig. 6.2 'e), may be associated to the oscillations 



of purity decay (inset of fig. |6.1| . 

Another example where the conditions of our derivation are not met is 
the following. A Bose-Einstein condensate in which the n atoms have 
two immiscible internal states |DD200| could be interpreted as a QM. 
The symmetry of the wave function reduces the dimension of the Hilbert 
space and causes high correlations among the couplings to the environ- 
ment. It is not surprising then that decoherence scales differently (as n^). 
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Figure 6.2: The color-coded abso' 



ute value of the real part of 



given in eq. (6.13 1 as a function of T and t' 



for the chaotic [(a) (b) (c)] and the integrable [(d) (e) (f)] spin 
chain environment (L = 8). (a),(d) show the autocorrelation 
function {i = j); (h), (e) and (c) (f) the cross correlation func- 
tion for couplings to the same spin or opposite spins in the ring 
respectively. 



Chapter 7 

Conclusions 



We now proceed to formulate the conclusior\s of this work. We start 
by particular conclusions for each chapter, and afterwards to the general 
conclusions and perspectives. 

In chapter |2] we analyzed the decoherence (measured by purity during 
all the thesis) of a single qubit immersed in an environment. For this 
purpose we introduced a model with unitary dynamics in the qubit plus 
environment space. Both the environment Hamiltonian and the interac- 
tion were chosen as random matrices from the classical ensembles. We 
solved purity decay in terms of correlation functions of the environment. 

We obtained linear and quadratic decay before and after the Heisenberg 
time respectively. In general, purity of a superposition of eigenstates 
decays faster than an eigenstate. The presence of an internal Hamiltonian 
in the qubit was observed to contribute to stability. For the GOE case, 
different initial conditions in the Bloch sphere yield different behaviors 
of purity and von Neumann entropy, even in the large environment limit. 

Monte Carlo simulations showed that exponentiating the analytical for- 
mulae, as explained in appendix |Cj allowed extending the analytical re- 
sult outside the scope of linear response formalism. 

In chapter |3] we used the spectator configuration (introduced in the in- 
troduction of the thesis), in which only a part (one qubit) of the central 
system (two qubits) is coupled to an environment. This constituted the 
first of three models analyzed. We solve it using extensively the sym- 
metries of the ensembles. We found some elements in common with the 
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single qubit case: the same qualitative dependence on time (linear and 
quadratic before and after the Heisenberg time, respectively; a depen- 
dence on an internal Hamiltonian, namely having a big energy splitting 
in the qubit increases stability; finally, a reasonable agreement with a 
heuristic exponentiation. With the second qubit we added a new ele- 
ment, entanglement. We observed that entanglement allowed taking the 
invariance properties of one qubit into the other. We also observed that 
entanglement induces a faster decoherence. For the TRI breaking case we 
found that, due to entanglement, self averaging is expected when using 
Bell states. 

The other two models (joint environment and separate environment), 
which possess a very different physical interpretation, were reduced to 
the spectator configuration using the independence of the couplings. Their 
behavior is similar, up to a scaling of the coupling constant. 

In chapter |4] we studied, for two qubits coupled to an environment, the 
relation between purity and concurrence as both evolve in time. Both the 
environment and the coupling are chosen from the classical ensembles. 
In order to analyze their relation we introduce the CP plane. We found 
that in the large dimension limit of the environment, concurrence and 
purity are restricted to a small region defined using unital channels. In 
fact, in the large dimension limit, for all parameters studied, the chan- 
nels induced by the random Hamiltonian approach rapidly the unitality 
condition. For large purities the allowed unital region is almost one di- 
mensional, which allows to write a one to one relation for concurrence 
and purity. In some cases the curve described by the system in the CP 
plane coincides with the Werner curve. The deviation from it, for a partic- 
ular level splitting in the coupled qubit, is given in an empirical formula 
which depends on the coupling strength and the size of the environment. 
This deviation goes rapidly to zero as the size of the environment and the 
coupling increase. Assuming that purity and concurrence are related via 
the Werner curve, it is possible to use our knowledge from the previous 
chapters to give a formula for concurrence decay. 

In chapter |5] we studied the same problem as in previous chapters (one 
and two qubit decoherence), but using particular dynamical models. The 
models consists of n spins interacting pairwise with individually ad- 
justable strengths, and kicked with a magnetic field. The dynamical sys- 
tem is flexible enough to allow different configurations and symmetries 
of the environment and coupling. Moreover, varying the parameters we 
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can have qualitatively different kinds of dynamics, namely chaotic, in- 
termediate and integrable. In the chaotic regime of the environment we 
obtained quantitative agreement with the RMT formulae obtained earlier 
for purity decay. For the integrable case we observed quadratic decay, 
sometimes with strong revivals (when the coupling was to an open end 
of a chain). The intermediate case displays oscillations superimposed to 
a linear decay. In the chaotic case we observe that the effective Heisen- 
berg time (determined, among other things, by the symmetry of both the 
environment and the coupling) plays an important roll in decoherence. 

Purity and concurrence behave in a qualitatively similar way for all con- 
figurations. The relation between the two, for sufficiently large envi- 
ronments, follows an analytic expression we get for Werner states quite 
closely, although the specific dynamics does not produce Werner states. 

Finally, in chapter |6] we calculated generic decoherence of a n-qubit quan- 
tum memory as represented by purity decay. For large purities it is glob- 
ally given in linear response approximation by a sum rule in terms of the 
purities of the individual qubits entangled with the remaining register as 
spectator. This sum rule depends crucially on the absence or rapid decay 
of correlations between the couplings of different qubits in the interaction 
picture. We prove that this is fulfilled by sufficiently chaotic environments 
or couplings but, using a spin chain model as an environment, we find 
that even if the latter is integrable correlations can be absent and the sum 
rule holds. While exceptions can be constructed we have a very general 
tool to reduce the decoherence of a QM of n qubits to the problem of a 
single qubit entangled with the rest of the QM. Furthermore the result is 
equally valid if we perform local operations on the qubits. 

Some possible extensions or generalizations of the thesis are now out- 
lined. 

• Several measures for entanglement of multi-qubit systems exist. De- 
pending on the particular problem to be treated, one or the other 
might be useful. Preliminary results |ste| indicate that for multi- 
qubit systems a relation, like the one found for two qubits, might 
be used to obtain formulae for entanglement decay. 

• In order to deepen our knowledge of the behavior of integrable and 
intermediate systems, one should study particular models. 

• As many interactions are indeed separable, a RMT formulation for 



84 



Chapter 7. Conclusions 



separable interactions of the qubit(s) with the environment might 
proof useful to understand the similarities and differences with this 
"more realistic" coupling. 

• The sum rule studied in chapter |6] also applies if the register is split 
into arbitrary sets of qubits. In fact, one and two qubit gates are 
known to be universal for quantum computation [NCOO|, so we lay 
the foundation for the computation of decoherence during the ex- 
ecution of a general algorithm. This extension only requires the 
knowledge of the decoherence of a pair suffering the gate operation 
while entangled with the rest of the register. As each gate is dif- 
ferent, we will have to take a step by step approach, which at each 
step will only involve one and two qubit decoherence. The linear 
response approximation will be sufficient due to the high fidelity 
requirements of quantum computation. 

Summarizing, we have developed a tool to study decoherence of qubits 
based on the spectator configuration and on random matrix theory. We 
have compared the results found with a kicked Ising model. We almost al- 
ways observed linear/ quadratic decay of purity before/after the Heisen- 
berg time of the environment and self averaging both with respect to the 
Hamiltonian and to the initial conditions. The exceptions are: (i) In the 
presence of an internal Hamiltonian, the central system eigenstates decay 
linearly even after Heisenberg time, (ii) If we have a time reversal invari- 
ant Hamiltonian, no self averaging with respect to the initial conditions 
is observed after the Heisenberg time, (iii) In the KI model, and presum- 
ably due to the separable structure of the Hamiltonian, no self averaging, 
before the Heisenberg time, was observed. Entanglement enhances de- 
coherence and may transport symmetry properties from one qubit to the 
other. When studying two qubits we also studied their internal entangle- 
ment, as quantified by concurrence. A one to one correspondence was 
often found, characterized by the Werner curve. Thus, formulae for en- 
tanglement decay are also available. For a larger number of qubits a sum 
rule can be applied in a wide range of circumstances to express purity 
decay of a group of particles in terms of purity decay of each particle in 
a spectator configuration. In all cases, heuristic exponentiation leads to 
good agreement beyond linear response scope. 



Appendix A 



One-qubit purity decay for 
pure and mixed states 



Here, we consider the one-qubit decoherence in a more general context. 
One-qubit decoherence as described in sec. |2] deals with a separable initial 
state Qo = pi'^ \^e){^e\, Pi = Thus Qq is a pure state of a single 

qubit coupled to an environment in the pure state \ipe)- We consider here 
arbitrary, finite dimensions of Hi and He, and we allow pi to be mixed. 
The latter poses some problems on the physical interpretation of purity 
as a measure for decoherence. However, these can be avoided by taking 
the point of view of the environment, as it becomes entangled with the 
central system. That means we consider the purity of the state of the en- 
vironment after tracing out the central-system's degrees of freedom; see 
eq. ( 3.10| and eq. ( |3.11 1. If we consider the entanglement with a specta- 



tor being responsible for the mixedness of pi, we can use the following 
results to describe purity decay in the spectator model. 

We will derive explicit expressions for the average purity {P{t)) as a func- 
tion of time, where the average is only over the random coupling. These 
expressions involve a few elemental spectral correlation functions, whose 



properties will be discussed in A.5 As detailed in sec. 2.2 we work with 
the echo operator M{t) in the linear response approximation eq. ( |2.9| . 
Note that I{t) is Hermitian, while, due to time-order inversion, J[t) is 
not. 

This section is slightly boring, but since this is a Ph.D. thesis, the details 
must be incorporated. My thanks an admiration to the devoted reader 
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who reached this section. 



A.l General calculation 



Let Qi and (?2 be two operators acting on the Hilbert space H' = He- 
We define the purity form p [ ■ ] as a function of pairs of such operators. 



P[Qi <^ Qi] = tre[tri(^i) tri(^2)] • 



(A.1) 



Since any Hnear operator acting on Ti' (S> Ti.' can be expanded in terms 
of separable operators of the form qi (g Q2, linearity implicitly defines the 
purity form for arbitrary linear operators. For arbitrary operators A,B 
acting both on H', the purity form has the following property: 



p[A (g) B] = p[A^ (g> B^]* = p[B (g) A] = p[B^ (g> A^] 



ti* 



(A.2) 



The purity defined in eq. ( |1.8| can be expressed in terms of the purity 
form as 

P{treQ) = p[Q0Q]. (A.3) 

Note that on the RHS of this equation, we first take the trace over the 
central system, i.e. we consider the purity of the state of the environment 
after tracing out the central degrees of freedom. With this little twist we 
can describe the one-qubit decoherence and the spectator model at the 
same time. 



To average purity, we use eq. (2.9 1 to compute Q^{t) Q^{i) in linear 
response approximation 



{P{t)) = {p[Q^{t)0Q''{t)])=p[{Q^{t) 

Keeping terms only up to second order in A: 

- iA Iqo ^Qo- QoI^ ^Qo + Qo^ IQo - Qq^ QqI 

- A^ 7^0 ^Qo + Qo}^ 0Qo + Qo^ }Qo + Qo^ Qo} 



(A.4) 



(A.5) 



A.l. General calculation 
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In the next step, we perform the ensemble average over the coupling. The 
odd terms in A vanish as purity (even without averaging) must be real. 
For the remaining terms, we use the properties in eq. ( A.2| and the fact 
that I{t) is Hermitian, to obtain 

{P{t)) = P{0) - A\Aj - Ai) (A.6) 

with P(0) = p[qo<^ Qo] and 

Aj = 4Rep[(/)^oO(?o] 

Ai = 2 (p[(%J) ^ Qo] - Rep[(I(?o (S) Iqo)] + ® Qol)]) ■ (A.7) 

In order to pull out the same double integral of Aj and Aj, we use the 
time ordering symbol T. It allows to write for the average purity as a 
function of time 

{P(t)) =P{0)-2\^ f dr f dr'ReAn (A.8) 

JO Jo 

with 

Aji = p[T{VV')qo ® Qo] - p[{V'qoV) ® qq] 

+ p[{Vqo V'qo)] - p[{V'qo ^ QoV)] , (A.9) 

where V and V' are short forms for the coupling matrices V{t) and V{t'), 
respectively. The arguments t and t' of the coupling matrices are inter- 
changed in the second and fourth term of Ajj. This does not change the 
value of the integral of course, but it facilitates to handle common terms 
in the following considerations. 



The integrand Aji in eq. ( A.9| consists of four terms. These terms will be 



considered, one after the other. We will always first average the argument 
of the purity form. Only then we perform the partial traces over subsys- 
tem Ti.-[, and therefore the final trace over the environment. The averaging 
is only over the coupling, and it is done by applying two simple rules, as 
described below. 

For the coupling Vi g in the product eigenbasis of Hq, we use either ran- 
dom Gaussian unitary (GUE) or orthogonal (GOE) matrices. Their sta- 
tistical properties are completely characterized by the following second 
moments (for notational ease we ignore the subscript i^e for a moment) 



{yijyki) = 5 a + agoe ^ji , 



(A.IO) 
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where Xgoe = 1 if V is taken from the GOE, while Xgoe = if V is taken 
from the GUE. In the interaction picture we then find: 

{Vij Vki) = { Su 5ju e-HE;-^^.) + xgo^ 3,^ Sj, e-(E,-E,) (r+r') ) _ ^^^^^ 

A.2 The GUE case 

We now procede to calculate one by one the contributions. 

p [T{VV') Qo ^o] (t)'- We first compute the average 

r {VV'qo ^Qo)= X] k;) {V{T)ij,kiV{T')ki,,„„){mn\Qo Qo 

ijklmn 

= E k7)e-^(^"-^")l"-"'l {ijlQo ® Qo, (A.12) 
ijkl 



where we have used the averaging rule, eq. ( A.ll[ l, as well as the fact that 



the time-ordering operator T requests to exchange t with t' whenever 
T < t'. The indices i, k and m denote basis states in Ti.^, while the indices 
/ and n denote basis states in He- We can rewrite that expression in a 
more compact form by employing the diagonal matrices Cj(t), defined 
in eq. ( |A.37| : 



r{VV'Qo(S)Qo) = (Ci(|T-T'|)®Ce(|T-T'|)) 

X {pi^\xpe){tpe\)^{pl(S)\xpe){%\)- (A.13) 

Now we apply the partial traces over subsystem Tii 

tr,{{VV')Qo) tri (?o = Ci(|t-t'|) Ce(|T - t'|) |t/;e) (t/'el . (A.14) 
The final trace over the environment yields 

p[T{VV')Qo(^Qo]it) = Ci(|t-t'|) CedT-r'l) . (A.15) 

pH'^'QoV) Qo] (t): We first compute the average 

ijklnmpq 

= L |2;)e-^(='>-^«)("-"')(W|()o|W)(/7| Qo, (A.16) 

ijkl 



A.2. The GUE case 
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where we have used eq. ( A.ll| . We apply the partial traces over subsys- 
tem Hi 

tri {{V'qoV)) tri qo = (^Ee-i(^'-^^)(---')(fc|pi|fc)^ (A.17) 

= Ci(t-t')EI;) e-i(^i-^'H^--') 

The final trace over the erivironmervt yields 

i' 

= Ci(T-T')Se(T-T'), (A.18) 

where Sj(t) is defined in eq. 



P[{Vqo (8) V'Qo)]it)- We first compute the average 

{V Qo (g) V' Qo) = I Qo \mn) Vmn,pq) {pq\ Qo 

ijkl mn pq 

= E e-'(=«-^'^) {kl\ Qo ® 1^/) Qo • (A.19) 

ijkl 

Then we apply the partial traces over subsystem Hi 

( tn (y(^o) tniv'Qo)) = (E^.e-^"*-"''^"-"') {k\pi\i) {i\pi\k) ) 

i' 

The final trace over the environment yields 

p[{VQo (S) V'qo)] = Si(t - t') Se(T - t') . (A.20) 

p[{V'qo0 QoV)]{t): We first compute the average 



{V' QqI^QoV) = E I '7) y!j,kl (^'1 Q0(^Q0\mn) Vnin,pq ^ 
ijkl mn pq 

= E 1^7) e-^(^'V-E«) (^/| ® ()o 1^/) (A.21) 
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We apply the partial traces over subsystem Hi 

{tri{VQo) triiQoV')) = ( E,,e-(^'-^^) (^-^'^ {k\pi\i) {i\pi\k) ) 

The final trace over the environment yields 

p[{VQo ^ QoV')] = Si(t - t') Ce(T - t') . (A.22) 

A.3 The GOE case 

In the GOE case, the average over the coupling yields, besides the pre- 



viously considered GUE-term an additional one; see eq. (A.lll. In the 
following we will redo the calculation of the previous subsection, but 
consider only that additional term. As a reminder, we add the subscript 
2 to the brackets which denote the ensemble average. 

p [T(yy') 2^0 ® ^o] (t)- We first compute the average 

r {V V')2 Qo(^QO= E i^i^)iiJd V{T')kl,mn)2 {mn\ Qo O Qo 
ij kl mn 

= El';> (27ko®(?o (A.23) 
'7 

= Qo(^Qo, (A.24) 



where we have used the averaging rule, eq. ( A.ll| |. Here the time-ordering 



operator has no effect. Now we apply the partial traces over subsystem 

tri ( {VV')2Qo) tri Qo = \tpe) {%\ • (A.25) 
The final trace over the environment yields 

p[{VV')2Qo^Qom = l. (A.26) 

^[(^'^0^)2 (8> Qo] (t)- We first compute the average 

{V' Qo V)2 '^Qo= E I'/) {^ii,ki {kl\Qo\mn) y,„n,pg)^ {pq\ ® 

ij kl mn pq 

= ^|z7)e-i(^'V-E«) {kl\Qo\ij) {kl\^Qo, (A.27) 

ijkl 



A.3. The GOE case 
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where we have used eq. ( A.ll| . We apply the partial traces over subsys- 
tem Hi 

trii{V'QoV)2) tn Qo = El7) e-(^;-^')(-+^') {l\ip,) {ip,\i) {l\tp,) {tp,\ . 

(A.28) 

The final trace over the environment yields 

i' 

= S;(-t-t'), (A.29) 



where S^(t) is defined in eq. 



p[{Vqo ® V'qo)2] (t)- We first compute the average 

{V Qo O V' Qo)2 = {^ij'ki {kl\ Qo \mn) Vm„,pq)2 {pq\ Qo 

ijkl nm pq 

= E e-^(^«-^") (^+^') {kl I Qo ® 1/7) {kl I Qo . (A.30) 

ijkl 

Then we apply the partial traces over subsystem TLi 

(tn {VQo)tr, {V'qo)), = (E,,ce-^('*-^')(^+^')(^|pi|0(^IPilO) 

X ^|y)e-i(E'-E.)(-+^')(/|^^)(^^|y)(/|^^)(^^|. (A.31) 

The final trace over the environment yields 

P[{Vqo^V'qo)2]=S[{t + t') S;(t + t') . (A.32) 

pH^'Qo ^ QoV)2] (t)- We first compute the average 



{V' Qo <S) Qo V)2 = E 1^'/) ylj,kl (^'1 ^0 ^0 1"^") ^mn,pq)^ {pq\ 
ijkl mn pq 

= E 1^'/) e-i(%-^«) (^+^') {kl \qo®Qo {kl I . (A.33) 

ijkl 

We apply the partial traces over subsystem Hi 

{tri{VQo) tn {qoV'))^= (E,,e-(='-^^*) {k\pi\i) {k\pi\i) ) 

xEli) e-(=i-^') (/IV^e) (V^el;) (/|. 
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The final trace over the environment yields 

p[{VQo ^ QoV')2] = S[{-T - t') . (A.34) 



A.4 The general solution 



If the coupling matrix is taken either from the Gaussian unitary (GUE) or 
orthogonal (GOE) ensemble, we find 



At 



[Ci(|t-t'|)-Si(t-t')] [Ce(|T-T'|)-Se(T-T') 



+ ^GOE [1 - S;(-T - t') +S[{t + t')S;(t + t') - Si(-T - t')] • (A.35) 

These expressions are derived from the previous two sections, for the 
GUE case in |A.2l for the GOE case in |A.3l The correlation functions 



Cx,Sx,C'^ and S[ with x S {l,e} are defined and discussed in A.5 



If the dimension of the Hilbert space of the environment goes to infinity, 
the corresponding correlation functions simplify, as discussed in A.5 We 
are then left with 

Aji= [Ci(|t-t'|)-Si(t-t')] C(|t-t'|)+xgoe[1-S;(-t-t')] • 

(A.36) 



A.5 Some particular correlation functions 



Averaging over the perturbation (i.e. the coupling) leads to expressions 
which may involve the diagonal matrix 

C;,(t) = ^|z) e-i(^*-^')^ (z| xGl,e. (A.37) 
ik 

Here, x denotes one of the two subsystems, either the qubit or the en- 
vironment. Evidently, the energies E^- are the eigenvalues of the corre- 
sponding Hamiltonian. In the derivations in this appendix, the expecta- 
tion value of Cx{t) with respect to the initial state Qo — Pi ^ Pe are of 
particular importance. These are denoted by 



Cx{t) =tr{Cx{T) px) xGl,e 



(A.38) 



A.5. Some particular correlation functions 
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In this work, = \xpe){tpe\ is always a pure state. 

We will also encounter another type of correlation function, which may 
be defined as follows 

(t) = tr, [ e-^«-^ " e^«-^ ] . (A.39) 

If the initial state px is pure, this quantity becomes the return probability. 
At T = it gives the purity of px- In the case of GOE averages, we also 
encounter the correlation function 

S; (t) = tr, [ e-i«-^ ^ Px e^«^ ^ pi ] . (A.40) 



In the main part of this paper, we focus on the limit, where the dimension 
of the environment(s) becomes infinite. In that case it makes sense to 
perform an additional spectral average over Hg and/ or Hgi. In the case of 
Ce{t) this yields 

(Ce(T)) = C(t) leC(T) = (Ce(T)) = 1 + 3{t/Th) - &f (t/Th) , (A.41) 

where b2^\t) is the two-point spectral form factor with time measured 
in units of the Heisenberg time th for the corresponding spectral ensem- 
ble. In the limit of large dimension Ne = dim('^^e) we find for the other 
correlation functions: 

Se(T) = tr (e-^»^^|i/;e)(t/^e|e^''^"|t/^e)(V'e|) (A.42) 

S;(t) = tr (e-i""^|i/;e)(i/^e|ei^^^|</^:)('/'el) • (A.43) 



For random pure states, as considered in the present work, both correla- 
tion functions are at most of order one at t = 0. For t > they drop 
very quickly and soon become of order N^^. This happens on the same 
time scale, where Ce(T) drops from values of the order Ng to values of 
the order one. In that sense we consider these correlation functions to 



contribute only 0(Ng ^) corrections to the result given in eq. ( A.36 i. 



A single qubit is a two level system. The most general pure initial state 
is given by \a) = |1) fli + |2) a2, with + |fl2p = 1- The most gen- 
eral mixed state is given by pi = Ai \a) (fl| + A2 \b) {b\, with Ai, A2 > 0, 
real, Ai + A2 = 1, and \a), \b) arbitrary pure states with {a\b) = 0. We 
will now investigate the behaviour of the different correlation functions 
ReCi (t). Si (t), and S[{t) as it depends on the initial state and the Hamil- 
tonian Hi = |l)£i(l| + |2)£2(2|. 
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Appendix A. One-qubit purity decay for pure and mixed states 



(a) Assume pi = \a) {a\. Then 

Ci(t) = Y^\ai\^e-'^^'^-^-^^ = l + cosAT + i(|fl2p - sin At, 

ik 

(A.44) 

holds, so 

Re Ci (t) = 1 + cos At A = E2 - Ei • (A.45) 

(b) For pi = Ai + A2 \b) {b\ we still have 

Re Ci (t) = 1 + cos At . (A.46) 

(c) Assume pi = \a) {a\. Then we find for Si(t) = s{\a);r) 

Si(t) = J2 |fld'l«/c|'e-i(^'-^*)^ ^ |^^|4 ^ |^^|4 ^ 2 |aip|a2|' cos At . 

ik 

(A.47) 

This expression only depends on the absolute values squared of the 
coefficients fli and fl2- Therefore we may parametrize them without 
loss of generality as ax = cos (p and fl2 = sin (p. We then find 

= g,p + i^-g(i>) cos At = |fli|'^ + |fl2|^ = 1-2 si"^ ■ 

(A.48) 

(d) For a general mixed state pi = Ai (a| + A2 \b) {b\ we find 

Si(t) = Afs(|fl);T) + A|s(|b);T) + 2AiA2 J^ai4bkbJ cos(Ei - Efc)T 

ft 

= A?s(|fl);T) + A2s(|b);T) + 

4AiA2|aip|biP(l -cosAt), (A.49) 

where we have used that {a\b) — aib\ + 02^2 — ^- Note that the co- 
efficients fli, fl2/ ^1/ ^2 niay be arranged into a square unitary matrix, 
and must therefore be of the following general form 

fa, ^i,/ ei^cos,^ ^V^^^. (A.50) 

V«2 02/ \— e '''^ sm ^ e cos (pj 

This shows that s( |a); t) = s(\b);r) = + (1 — ^(^i) cos At, and that 

Si {r) = {X\ + \l)[g^ + {\- g^) cos At ] 

+ 2AiA2(l-^<^) (1-cosAt) . (A.51) 



A.5. Some particular correlation functions 
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Since Ai + A2 = 1 it is convenient to set Ai = cos'^ 9 and A2 = sin^ 9 
such that we may write A^ + A^ = ge and obtain 

Si{T) = l-g0-g^ + 2geg^ + i2ge-l)il-g^) cos At. (A.52) 
(e) Assume again that pi = \a) {a\. Then we find for S[{t) = s'(|fl);T) 

= + 2 |fli p |fl2 p cos( At + Irj) 
= 8<p + {'^- S<p) cos(At + lrj) t] = arg(fli) - arg(fl2) • (A.53) 



However, as discussed in sec. |2.5} a natural symmetry around the y 
axis (in the Bloch sphere picture) appears for A = 0. In that case, 
using 7 as defined in that section, one can prove that 

1 - Si (t) I^^^ = (1 - [1 - cos(2//)] = sin2 7 (A.54) 

using elementary geometric and trigonometric considerations. 

(f) For a general mixed state pi = Ai \a) (fl| + A2 \b) {b\ we find 

S[{t) = A?s'(|fl);T)+A2s'(|&);T) + 2AiA2Re J^a^albPie-'^^'-^^^^ . 

ik 

(A.55) 

It follows from eq. ( A.50| that arg(&i) — arg(&2) = arg(fli) — arg(fl2) — 
71, such that the equality s'(|fl);T) = s'(|&);t) holds, just as in the 
case of Si(t). Therefore we may write 

S[{r) =gg[g^ + {l- g^) cos(At + Irj) ] (A.56) 

+ 2 A1A2 (1 - Re( 1 + ei^^^+^-z-^) ) 

= ^-ge-g,p -^2g6g<p + {2ge - - g<p) cos(At + 2//) . 

(A.57) 

Note that the only difference to S(t) in case (d) is the additional 
phase Irj in the argument of the cosine function. Using the same 
angle 7 defined with eq. ( 3.26| l, in complete analogy with eq. ( 2.23| , 
and using eq. ( A.54| we obtain 

1 - Si(T)|^^Q = l-ge + {2ge - l)sin2 7. (A.58) 



Appendix B 



Implementing the evolution of 
the KI models 



In this appendix we describe how to implement numerically Hamilto- 
nians describing a set of qubits with pairwise variable Ising interaction, 
and kicked with a periodic, and site dependent, magnetic field. We have 
used this systems in chapters |5] and |6| the reason being (i) the flexibility 
of the dynamical model to describe various physical situations, (ii) it has 
various dynamical regimes, and (iii) it is has an efficient implementation 
making it useful for numerical studies. 

Recall eq. ( |5.5| l, the Hamiltonian of the system is 

;>/c=l ;=1 

which gives rise to the Floquet operator Uki = LTising^^kick with 




However notice that 

falsing = Yl ^^P(~^jhk'^'^k) ' U'kick = n^'^P (~i^r (^-3) 

j>k=l ;=1 

Thus the evolution operator is decomposed trivially into single and two 
qubit operators that can be efficiently implemented in a computer. Notice 
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Appendix B. Implementing the evolution of the KI models 



that we do not need to explicitly diagonalize or even store the complete 
evolution operator thus allowing the implementation of twice as much 
spins than with brute force diagonalization. 

The state can be stored in a complex array with indices ranging from to 
2^ — 1. The state (expressed in the computational basis) 

1 

1^) = E ^iL-iiL-2-ioh-lk-2---iQ), (B.4) 
iL-i,iL-2,- ,io=0 

can be stored in a complex array of dimension 2^, if in the position ij2^ 
we store the value */t_i/Y_2 'o- ^ Fortran 90 routine implementing the 
Ising interaction with the state stored in this form is the following. 



subroutine Action_Ising_Interaction(psi , 1 ,iii, J) 
coniplex(kind(ldO)) ,intent(inout) :: psi(0:) 
real(kind(ldO)) ,iiitent(iii) :: J 

integer, intent (in) : : l,iii 

integer : : i 

do i=0,size(psi)-l 

If (btest (i , 1) . eqv.btest (i ,m) ) then 
psi(i)=exp(-(OdO,ldO)*J)*psi(i) 

else 

psi(i)=exp((OdO,ldO)*J)*psi(i) 
end If 
end do 

end subroutine Action_Ising_Interaction 



A routine to implement the magnetic kick is now presented. 



subroutine Application_Of _Kick(psi ,b , 1) 



complex (kind (IdO) ) , intent (inout) 
real (kind (IdO)) .intent (in) 
complex (kind (IdO) ) 
complex (kind ( IdO) ) 
integer , intent (in) 
integer 
integer 
if (sum(b**2) . eq. OdO) return 
Matrix_For_Kick=Matrix_For_Magnetic_Kick(b) 



psi(0:) 
b(:) 

Matrix_For_Kick(2 , 2) 
z(2) 
1 

j _lef t , j _right , lef t_digit 
UpperNumber , LowerNumber 
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do j_left=0,size(psi)/2**(l+l)-l 
lef t_digit=ishf t ( j _lef t , 1+1) 
do j_right=0, 2**1-1 

LowerNumber=lef t_digit+j _right 

UpperNuinber=IBset (LowerNumber , 1) 

z(l)=psi(LowerNumber) ; z(2)=psi(UpperNumber) 

z=MatMul(Matrix_For_Kick,z) 

psi (LowerNumber) =z(l) ; psi(UpperNimiber)=z(2) 
end do 
end do 

end subroutine Application_Of _Kick 

The previous routine uses the following function. 

function Matrix_For_Magnetic_Kick(b) 

complex (kind ( IdO) ) :: Matrix_For_Magnetic_Kick(2,2) 

real (kind(ldO)) .intent (in) :: b(3) 

real (kind (IdO) ) :: n(3) ,phi,cosphis2,sinphis2 

phi=norma(b) ; n=b/phi 
Matrix_For_Magnetic_Kick(l , 1)=& 

cos(phi)-(OdO,ldO)*n(l)*sin(phi) 
Matrix_For_Magnetic_Kick(2 , 2) =& 

conjg(Matrix_For_Magnetic_Kick(l , 1) ) 
Matrix_For_Magnetic_Kick(l,2)=& 

-(n(3)+(0d0,ld0)*n(2))*sin(phi) 
Matrix_For_Magnetic_Kick(2, 1)=& 

-conjg(Matrix_For_Magnetic_Kick(l ,2) ) 
end function Matrix_For_Magnetic_Kick 



Appendix C 

The exponentiation 



The linear response formulae for purity derived throughout the thesis 
are expected to work for high purities or, equivalently, when the first 
two terms in the series eq. ( |2.9| approximate the echo operator well; it 
is of obvious interest to extend them to cover a larger range. The exact 
solution for fidelity decay has, in some cases, been possible with some 
effort, using super-symmetry techniques. This has been partly due to the 
simple structure of fidelity, but trying to use this approach for a more 
complicated object such as purity seems to be out of reach for the time 
being. Exponentiating the formulae obtained from the linear response 
formalism has proved to be in good agreement with the exact super- 
symmetric and numerical results for fidelity, if the perturbation is not too 
big. The exponentiation of the linear response formulae for purity can be 
compared with Monte Carlo simulations in order to prove its validity. We 
wish to explain the details required to implement this procedure in this 
appendix. 

Given a linear response formula Plr(0 t^o^ which P(0) = 1], and an 
expected asymptotic value for infinite time Poo, the exponentiation reads 

Pelr(0 =-Poo + (l-Poo)exp 

This particular form guaranties that 

Pelr(0 ~ PLR{t) (C.2) 

for short times, and that 

lim PELR(f) = Poo. (C.3) 



1 - Plr(0 



(C.l) 
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Appendix C. The exponentiation 



The particular value of Poo will depend on the physical situation; in our 
case it will depend on the configuration and on the initial condition. We 
now turn to an estimation of this quantity for the spectator configuration, 
under some very general assumptions. 

Let us write the initial condition in the qubits as 

\xp12i0)) = cos0|OiO2) +sin0|iii2) (C.4) 

for some rotated qubits |0/), We assume that for long enough times, 
decoherence is described for the uncoupled qubit by the following chan- 
nel: 

1 1 

£td[p] = 2^4 i'^^P^^ + '^yP'^v + '^^P'^^ + P) ■ (^-5) 

This is the depolarizing channel with full strength IINCOOI . This channel 
is unital. We now evaluate S^d (8) t[\tpi2{9)){^pi2{0)\]. Note that Std 
1[|00)(00|] = l^|0)(0|/2 and ftd«'l[|ll)(ll|] = lO|l)(l|/2. For the 
cross terms we obtain: 

£....1||00)(11|||= l'°><°^l'l'°><"l;l°°><"l + l°°><"l =0, (C.6) 

and also £^td (00| |] =0- Using the linearity of the channel and the 

results found, we see that 



In the end, we estimate the asymptotic value of purity as 

cos^ 6 + sin^ 9 _ ge 

^00 = ^ - y, (^-.Oj 

where ge was defined in eq. ( 2.17| |. 

For both the joint and the separate environment configuration we use the 
value 

Pco = P{Sd(^Sd[\tpi2{e)){ipi2{e)\]) = \. (C.9) 

Good agreement is found with Monte Carlo simulations for moderate and 
strong couplings. Obtaining the value of Poo for weak coupling should be 
possible using perturbation theory on the eigenvectors. This calculation 
is still missing. 
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Appendix D 



Entanglement 



Here we want to discuss or comment on some technical aspects that do 
not need to be in the introduction. We prefer to include them in an 
appendix to ease the reading of the main part of the text. 

Recall from eq. that a pure state of a bipartite system is separable if it 
is the tensor product of states in the individual systems. The generaliza- 
tion for multipartite systems is straightforward. Let us have n particles, 
each in a Hilbert space Tii, f = 1, ■ ■ ■ , n. A state \-ip) is said to be separable 
if it can be written as 

\^)=^\^r)^ m^^i), (D.l) 
1=1 

compare with eq. Otherwise it is called entangled. A mixed, bipar- 
tite system in state p is said to be separable if it can be written as a convex 
combination of separable states. That is, if 

P = T.P^pf^®P?^ (D.2) 

i 

where p,- > and p\^^ ip\^^) are density matrices of system A (B). Though 
this definition is simpler from a conceptual point of view, often a more 
comfortable formulation in terms of pure states if obtained expressing the 
density matrices in terms of projectors. Thus a bipartite mixed state p is 
said to be separable if it can be written as 
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Appendix D. Entanglement 



with \xpA) ^ Ti-A, IV'b) £ 'Hb and p, real and positive. For the multipartite 
case one should require, instead of eq. ( D.2| , that 



(k) 

with p- being a density matrix of space Hk and the piS probabilities {i.e. 
still a convex combination of separable states). 



D.l Quantifying two qubit entanglement 

Here we want to present the entanglement of formation, and its relation 
to concurrence. 

Let us first define the entanglement of formation. The entanglement of 
formation Ep measures the number of Bell pairs needed to create an en- 
semble of pure states representing the state to be studied IIHW97I . We 
define Ep by stages. The entanglement of formation for a pure state is 
the von Neumann entropy of the reduced density matrix. Given \xp) with 

Schmidt coefficients Ai and A2 = yjl — \\, Ef{\tp){xp\) = — AilogjAi — 

A2 log2 A;^ For mixed states we have 

£f(p) = min J2Pi^Fi\^i){^t\) such that p = Y2Pi\^i)i^i\ (^-5) 

{piM'i}} i i 

and the p,s are positive. This represents the minimum amount of entan- 
glement needed to create an ensemble that reproduces p. All the require- 
ments discussed above are fulfilled. 

Concurrence C is defined via the entanglement formation: 



£^(p) = /, ^1±V1_C(pP j . ^(^^ ^ -Xl0g2(x)-(1-X)l0g2(l-X). 

(D.6) 

This weired definition is inspired in the fact that a simple and close for- 
mula for C{p) was obtained. 



Notice the close relation with the classical entropy. 



D.2. Some generalizations 
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D.2 Some generalizations 

Two possible generalizations for the entanglement measures would be de- 
sirable: (i) Allow qudits instead of qubits (i.e. having more that 2 levels). 
Furthermore one could like to quantify the entanglement in systems with 
infinite or a continuous number of levels, (ii) Allow the Hilbert space of 
the system to have a more complicated partition. Namely let H = ^^^^Hk 
with n being bigger than 2; for example if one wants to examine the en- 
tanglement between three particles, n = 3. The cases in which the state 
representing the system is pure are again much simpler that when the 
state is mixed. 

For both cases the concept of entanglement is easy to define. A relevant 
review of entanglement measures for qudits is [P V07J . Some measures of 
entanglement, for the multipartite case, are discussed in IIMCKB 05I. As 
the reader can find, this generalizations are not simple, unique or easy to 
calculate. 



Appendix E 

RMT: various aspects 



In this appendix we wish to discuss some technicalities about the random 
matrix ensembles used in this thesis. It includes a small discussion about 
the ensembles with /without time reversal symmetry, the construction of 
the ensembles, and concepts as the Heisenberg time, level density, and 
form factor. 

We start by discussing the difference between the unitary ensembles (GUE 
and CUE) and the orthogonal ensembles (GOE and COE). We first de- 
scribe the GOE. This is an ensemble of Hermitian matrices with the spe- 
cial property of being real. As we know, any Hamiltonian describing a 
system with an anti-unitary symmetry (as time reversal) can be repre- 
sented as a real matrix without diagonalizing it. The COE corresponds to 
the unitary version of this ensemble, for each H in the GOE, a correspond- 
ing exp(zH) is found in the COE. On the other hand, when the physical 
system has no underlying anti-unitary symmetry, it must be described by 
arbitrary Hermitian matrices. The GUE is the appropriate ensemble as it 
has no anti-unitary symmetry. The set corresponding to this ensemble is 
the set of all Hermitian matrices. The CUE is the unitary version of this 
ensemble. The most common example of this distinction involves a parti- 
cle in a chaotic billiard. With /without a magnetic field the level statistics 
resemble those of the GUE/ GOE. Numerical l|BGdS95l and experimental 
fSA0095t .DDG+031 confirmation has been reported. 

We now turn our attention to the construction of the ensembles used in 
this thesis. The neatest way to perform that construction is using the 
arguments of least information (see llBal68l Dys62 1 or, for a more modern 
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Appendix E. RMT: various aspects 



view, IIMel95l section 2). We wish to recall that an ensemble is a set and 
a probability distribution, which can be regarded as a measure. The set 
shall be no problem to define, but the measure is more delicate. 

We start by defining the entropy of a distribution of matrices (objects) A. 
Let P(A) be a probability distribution on a set matrices A (or of anything). 
The information of the ensemble is measured via the information entropy 
HSha48l : 

J[P(A)] = - ^ P(A)logP(A)dA (E.l) 

Notice that we need an a -prior measure dA. This measure, in all the cases 
considered here, is an invariant (under a suitable operation) measure. If 
there is a group structure underling (like for the GOE and GUE with 
matrix addition and the CUE with matrix multiplication) this measure 
coincides with the Haar measure. Otherwise (for the COE), it is slightly 
more complicated. For details see IICar35l lMeh91 1. The information en- 
tropy, for a given can be interpreted in a number of ways [Sha48J. 
For example, it measures the amount of information in a given member 
of the ensemble. For bigger entropy, the information in each member is 
smaller {i.e. its members are closer to the anonymity). 

Let us start by constructing the CUE, which stands for circular unitary en- 
semble. This is an ensemble of unitary matrices of given dimension. The 
distribution that we use for the ensemble is simply the homogeneous one, 
P(A) = constant. Intuitively it is clear that, for this set, it is the ensemble 



with greatest information entropy, see eq. (E.l I. We can do an analogous 
treatment for orthogonal matrices to obtain the circular orthogonal en- 
semble (COE), using a uniform "probability" for all orthogonal matrices. 

For Hermitian operators, as compared to the unitary ones, the situation is 
more complicated. Think of Hermitian 1x1 operators (real numbers). Its 
Haar measure is simply the usual Lebesgue measure. The whole set has 
an unbounded measure hence it cannot be normalized; working naively 
with oc 1 results in diverging integrals. Moreover, for the real line 

(and using the usual measure) there does not exists such a thing as a uni- 
form distribution. Similar problems result for larger matrices due to the 
non-compactness of the set under consideration. However one can create 
a least information ensemble if one is willing to postulate an additional 
condition. We shall require that the average of tr (the "strength" of the 
operator) is some fixed number. Its particular value is irrelevant and will 
only impose some normalization condition in the ensemble. Notice that 
this condition is basis independent. One can solve the problem of max- 
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imizing the entropy constrained to the condition above, using Lagrange 
multipHers. The solution is 

P(H) oce-^''"^' (E.2) 
with A a real (ensemble dependent) number. 

In all the cases (COE, CUE, GOE, and GUE) we can see that the distri- 
bution is invariant under unitary (for the CUE and GUE) or orthogo- 
nal (for the COE and GOE) transformations, e. g. for the GUE, P{H) = 
P{U^^HU) for any unitary operator U, as the only condition we imposed 
is basis independent. The invariance properties of the ensemble will be 
used heavily to simplify the problem under consideration and group to- 
gether equivalent cases during the thesis. 

A practical way to realize the Gaussian ensembles (G*E) is to choose the 
matrix elements as Gaussian random variables with a given standard 
deviation a. Consider the independent elements of a Hermitian matrix 
H. We have that trH^ = E; + 2£,<y H,yHy; and then the distribution 
eq. ( |E.2| can be factorized as 

N N 

P(H)oc Yl e-2AH,«;'ne-^«i. (E.3) 

i<j=l ;=1 

For the GOE (of fixed dimension N) we choose: the diagonal elements 
as real variables with = 1 (forget about the global, unimportant nor- 
malization constant); for the elements above the diagonal, real variables 
with (7^ = 2; and the others (below the diagonal) are chosen to make the 
matrix symmetric. The matrix elements average like 

Wj = 0, VijWi = ^n^jk + SikSjj . (E.4) 

Of course the matrix elements can be multiplied by a convenient (possibly 
N dependent) factor to achieve certain particular condition. Under sim- 
ilar considerations one can show that a typical GUE element is obtained 
as follows. Choose the diagonal elements from a real Gaussian ensemble 
with cr = 1. The upper triangle of the matrix are complex numbers; the 
real and imaginary parts are chosen independently from a real Gaussian 
ensemble with a = 1/ \/2. The lower elements are chosen to make the 
matrix Hermitian. Its matrix elements average as 

Wj = 0, W^i = &ii&jk- (E.5) 
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Figure E.l: Visualization of the Zeno time (above) and the Heisen- 
berg time (below). Eigenvalues E are plotted in the unit circle 
as dots as time passes, using exp{—iEt/h). An arc between 2 
particular adjacent eigenvalues is shown in red. The Zeno time 
Tz is close to the first time chosen all eigenvalues cover the unit 
circle. Heisenberg time th is the average time when adjacent 
eigenvalues cross in the unit circle. 100 Poissonian eigenvalues 
were used for these plots. 



At this point we want to introduce the concept of level density, sometimes 
also called density of states. This is the number of levels per unit interval. 
The level density is a particular feature of each system, even of chaotic 
ones. Random matrices have their own level density. The circular ensem- 
bles have a constant level density (within the interval [0, 27r)) whereas the 
Gaussian ensembles have an ellipsoidal level density (often referred to as 
circular). In the large dimension limit it is given by 



for the normalization previously proposed, see fig. E.2 Other concept 
closely related to the level density is the Heisenberg time Th. This is the 
typical time in which two neighboring levels cross in the unit circle [i.e. 
when (E,+i — Ei)t/fi In). Of course if the level density is not uniform, 
the Heisenberg time is also not uniform. Sometimes it is convenient to 
have a uniform Heisenberg time. This amounts to have constant density 
of states throughout the spectrum. If one uses the circular ensembles this 
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is automatically granted. Otherwise, the constant density can be achieved 
using the unfolding procedure. Normally it requires adjusting an average 
density but, since for the GUE and GOE we know the average density, we 
can calculate its cumulative function |GMGW98I , and use that to unfold 
the spectrum. Let Ec be the smallest expected energy (eigenvalue) of 
your spectrum. The urvfolded energy cpi corresponding to the original 
eigenvalue E, is 



N 
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sm ^ e,- + e,- 



Ec' 



(E.7) 



Another important time scale is the Zeno time Tz- It is the time the eigen- 
values are far from the origin, but are still much smaller than 2n. Both 



concepts are sketched in fig. E.l 



An important characteristic of the RMT spectra are its correlations. We 
can use the form factor 



K2{t) :-- 



1 

N 



N 
7=1 



(E.8) 



to quantify this correlation. In the large dimension limit we have that 



K2it) = 4f^\t) = 1 + 3 



t 



t 

Th 



(E.9) 



The index /3 characterizes the ensemble, /5 = 1 for the GOE (and COE) 
and j6 = 2 for the GUE (and CUE). The functions hi^"^ are given by 




and 



,(2), 



ln(2|i| + l) if|i|<l. 



if \t\ > 1 

if |f| < 1, 
if Ifl > 1 ' 



(E.IO) 



(E.ll) 



see fig. E.2 



A characteristic features of the GOE and GUE spectra is the hole in the 
correlation function (for a random spectra, K2{t) = 1 + S{t/TH))- This 
hole has been shown to induce a certain stability in fidelity decay. In this 
work it will also display some enhancement of stability though not as 
spectacular as for fidelity decay. 
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Figure E.2: We see here the level density of an ensemble of 800 ele- 
ments, of the GOE (top) and GUE (bottom) with N = 50. We 
show, from left to right, the level density eq. (|E.6|, the unfolded 



level density (see eq. (E.7|) and the form factor K2 Eqs. |E.10 and 



E.ll 



Appendix F 



On the numerics of RMT and 
quantum information 

We now describe some tricks to implement some aspects of RMT. We also 
include some of the routines useful for quantum information simulation. 
All routines are for Fortran 90. pi = n and 1 = i. 



F.l Random matrices 

To build the random matrices, as noted in Appendix |E| one can use ran- 
dom Gaussian variables. Most low level languages provide homogeneous 
distributions in the unit interval [0, 1 ) . A transformation from this dis- 
tribution (characterizing variables u and v), to a complex Gaussian one 
(characterizing variable z) of width a and centered at xq is provided by 

z = (je^^'"" 7-2 log M + xo (F. 1 ) 

Its implementation is straightforward, however if w = (which occurs nu- 
merically with very small, but finite probability) the routine could crash. 
If instead of selecting u G [0,1) we select 1 — u G [0,1), the problem is 
fixed. The routine is called RandomGaussian. Its first argument is a and 
the second one xq. 

In view of the considerations of the Appendix |E| to build a GUE matrix, 
one has to choose the diagonal elements real and with = \ whereas for 
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the off diagonal elements tr^ = 1/2. The routine to build the GUE matrix 
is: 

subroutine gener at e_GUE_inatrix (matrix) 
implicit none 

integer :: nsize,jl,j2 

complex (kind ( IdO) ) , intent (out) :: matrix(:,:) 
nsize=size (matrix ( 1 , : ) ) 
do jl=l,nsize 

matrix (j 1 , j l)=real(RandomGaussian(ldO,OdO) ) 

do j2=jl+l,nsize 

matrix (j 1 , j2)=RandomGaussian(l/sqrt (2d0) , OdO) 
matrix(j2, j l)=conjg(matrix(jl , j2) ) 

end do 
end do 

end subroutine generate_GUE_matrix 



A routine to generate a GOE matrix is 



subroutine generate_GOE_matrix (matrix) 
implicit none 

integer :: nsize,jl,j2 

real (kind ( IdO) ) , intent (out) :: matrix (:,:) 
nsize=size (matrix (1, :)) 
do jl=l,nsize 

matrix (j 1 , j !)=& 

real(RandomGaussian(sqrt(2dO) ,OdO) ,kind(ldO)) 
do j2=jl+l,nsize 

matrix (j 1 , j2)=real (RandomGaussianO ,kind(ldO) ) 
matrix(j2, j l)=matrix(j 1 , j2) 
end do 
end do 

end subroutine generate_GDE_matrix 



The normalization is such that eq. (E.4| holds. To obtain the spectrum 



of such matrices one must diagonalize numerically a dense matrix. The 
LAPACK90 package provides this routines conveniently packed for Fortran 
90. 



Most RMT results refer to the center of the spectrum, where the density 
of states is constant. To have better agreement of the numerics with the 
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theory one can unfold the spectrum, i.e. force an approximate constant 
level density. Given a density p{E), the function rj{e) = J^^^^^^dEp{E) 
maps £ to an "unfolded" energy e. Since 

^ yi-£2 = ^ ( e _ e2 + 71 - sin-2 (F.2) 

(which would assume an spectrum with ranging roughly from ±1), the 
following routine gives an unfolded spectrum, ranging from to N. 

function approxcumulative (Energy , size_spectrum) 

real (kind ( IdO) ) :: approxcumulative 

real (kind ( IdO) ) , intent(in) :: energy 

integer, intent (in) : : size_spectrum 

real (kind ( IdO) ) :: rescaled_energy 

rescaled_energy=Energy*sqrt (pi) /size_spectrum 

if (rescaled_energy <= -IDO) approxcumulative=OdO 

if (rescaled_energy >= IDO) approxcumulative=pi 

if ((-ldO<rescaled_energy) .and. (rescaled_energy<ldO)) & 

approxcumulative= (asin(rescaled_energy) + & 
rescaled_energy*sqrt (l-rescaled_energy**2)+pi/2) 
approxcumulative=(approxcumulative-pi/2) *size_spectrum/pi 

end function approxcumulative 

We assume in this function that the spectrum ranges from ±N/ \fn. 
The unfolding procedure can also work with the same routine presented 
above as the level density is again a semicircle. 



F.2 Quantum information 

During the numerical implementation of the ideas exposed in the the- 
sis, we have to do operations that involve tensor product of states [as in 
eq. (|Z4)], like 

\^) = \^a) ® \tpB) (F.3) 

with \ipA) G Ha and Ii/'b) G Hb- We shall assume that both Ha and Hb 
are spaces of one or more qubits. This allows to take advantage of the 
bitwise operations available in all programing languages. 

In order to explain the implementation of tensor products its better to 
use an example. We restrict our selfs to tensor products of qubit spaces. 
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Consider an space composed of L qubits, numbered from to L — 1: 
Ti. = Ti-i-i (8) • • • (8) Ho- Now consider the two subspaces Ha = ^aa-i ® 
■ • • (8) Hao of dimension 2"^ and Hb = 'Hhg_^ (8 • • • (8 Hhg of dimension 
2^ such that Ti. = Ha ®Hb- Assume that L = 7, Ha is the subspace 
of qubits 0, 3 and 5 and thus Hb is the subspace of qubits 1, 2, A, 6. 
One can label unambiguously Ha with a number tia = J^i^."'. In our 
case, ria = 2*^ + 2^ + 2^ =41. Analogously one can label Hb with ri}, = 
2^ + 2^ + 2"^ + 2^ = 86 (notice that ria + m = 2^-1 and thus are not 
independent). Using tensor notation, one can associate an single index fi 
of H with two indices, one in Ha and one in Hb- If one uses the notation 
explained in Appendix |B] one must combine the binary bits of index ji 
with the bits of the numbers representing each space (say /a in Ha and 
z'b in Hb) ■ As an example assume that f/ = 55. One obtains = 5 and 
= 7, as shows the following table. 



qubit number 


6 5 4 3 2 1 


f/ = 55 


110 111 


n„ =41 
lA = 5 


10 10 1 
10 1 


ni, = 86 
/b = 7 


10 10 110 
1 11 



A routine implementing the procedure is 

subroutine exdig(number_in,L,nlout ,n2out ,nwhich) 
implicit none 

integer, intent (in) : : L,nwhich,number_in 
integer, intent (out) : : nlout,n2out 
integer :: j,nuiiiin, L1,L2 

nlout=0 
n2out=0 

num i n=nunib e r _ i n 
Ll=bits_on_one (nwhich) 
L2=L-L1 
do j=0,L-l 

if (btest (nwhich, j ) ) then 

if (btest (numin, 0) ) nlout=ibset (nlout , 0) 
nlout=ishf tc (nlout , -1 ,L1) 
else 

if (btest (numin, 0) ) n2out=ibset (n2out , 0) 
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n2out=ishf tc (n2out , -1 , L2) 
endif 

inimin=ishf tc (numin , -1 , L) 
end do 
return 
end subroutine exdig 

Here number_in = ^, L = L, nlout = Ia , n2out = Zg, and nwhich = ria- A 
routme to calculate tensor products of real matrices, that takes advantage 
of the previous routine is the following. 

subrout ine TensorProduct _real_rout ine (H_a , H_b , nwhich , H_out ) 
integer, intent (in) : : nwhich 

real (kind ( IdO) ) , intent (in) :: H_a(0: ,0: ) ,H_b(0: ,0: ) 
real (kind ( IdO) ) , intent (out) : : H_out(0:,0:) 
integer :: total_size, jl, j2, & 

jla, jib, j2a, j2b, qubits 
total_size=size (H_a, l)*size(H_b, 1) 
qubits=IntegerLogBase2(total_size) 
do jl=0,total_size-l 
do j2=0,total_size-l 

call exdig (j 1 , qubits , j la, j lb , nwhich) 
call exdig (j 2 , qubits , j 2a , j 2b , nwhich) 
H_out(jl, j2)=H_a(jla, j2a)*H_b(jlb, j2b) 
enddo 
enddo 

The function IntegerLogBase2 calculates the integer logarithm in base 2, 
of its argument (it must be a power of 2). Again nwhich = n^, H_A (H^) is 
the operator acting in Ha CHa)- 

The last routine that we wish to document is one for performing partial 
traces. This routine calculates the partial trace over Hb of a pure state 
p = trg \tp){^\- The input parameters are the state \ip) (statin) and 
(nwhich). The output parameter rho is p. 

subrout ine safe_PartialTrace(statin, rho , nwhi ch) 
implicit none 

complex (kind (IdO) ) :: statin(:) 

complex (kind (IdO)) :: rho(0:,0:) 
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integer, intent (in) 
integer 



: : nwhich 
: : & 

err (3), n_final, jl, j2,qubits 



complex (kind ( IdO) ) , pointer :: & 



staout (:,:), stl(:), st2(:) 



qubits=IntegerLogBase2 (size (statin) ) 
n_f inal=bits_on_one (nwhich) 

allocate (staout (0 : 2** (qubits-n_f inal) -1,0: 2**n_f inal-1) ,& 

stat=err(2)) 
allocate (stl(0:2**(qubits-n_f inal)-l) ,& 

st2(0:2**(qubits-n_f inal)-l) , stat=err(3)) 
if (any (err /= 0)) then 

print* , "error en PartialTrace" 
stop 
end if 

call unmsta(statin, staout , nwhich, qubits) 
do j 1=0, 2**n_f inal-1 

do j 2=0, 2**n_f inal-1 

rho(j 1 ,j2)=dot .product (staout ( : , j2) ,staout(: , jl)) 

enddo 
enddo 

deallocate (staout, stl, st2, stat=err(l)) 
return 

contains 

subroutine unmsta(statin, staout , nwhich, qubits) 
implicit none 

integer, intent(in) : : nwhich, qubits 

complex (kind (IdO) ) , intent(in) :: statin(:) 
complex (kind ( IdO) ) , intent(out) :: staout(:,:) 
integer :: j, ncol, nrow 

staout=OdO 

do j=0,size( statin)-l 

call exdig(j , qubits, ncol, nrow, nwhich) 
staout (nrow+1 ,ncol+l)=statin(j+l) 
enddo 
return 
end subroutine unmsta 
end subroutine saf e_PartialTrace 

The function bits on one counts the number of bits of its argument, set 
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to 1. The subroutine unmsta reorganizes the input state as a matrix, to be 
able to express the partial trace as a dot product of the vectors composing 
this matrix. All the programs, and routines are available upon explicit 
request to the author. 



Appendix G 

Two minor technicalities 



G.l The Form Factor 



We define the spectral form factor as 

1 N 

m = ^Lf^'- (G.l) 

Some authors rescale the energy (or, equivalently the time). The factor in 
front is used to obtain an asymptotic average value of one for most spec- 
tra. This quantity has been studied extensively in the context of quantum 
chaos LGMGW98.JPAJ 

We want to evaluate a double integral of the two-level form factor for the 
GUE case: 

Let us suppose first that < t < th, and thus replace &2 by a linear 
expression: 

Bm = l.rl'.r',r{^) ,G.3, 

l-I^) ,G.4, 

2 - W„ 



f dT r dr' 
Jo Jo 
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If t > th the expression must be strictly linear since b2{t > th) = 0, and 
with slope given by J^^ h\[t/ th) = th/2. Furthermore, since B2{t) must 
be continuous, we also have the y intersect of the linear function. This 
leave us with the result 



2 6^ if < f < Th 
Bm={ 1 ■ (G.6) 



2 6 



G.2 A small proof of the Born expansion 

Consider the Hamiltonian = Hq + AW. The definition of the state ket 
in the interaction picture is IISak94| 

\^p{t))i = Ul{t)U,{t)\^{Q)) (G.7) 

where U\{t) = exp{itH\). Consider the equation of motion of the ket in 
the interaction picture IISak94| : 

r^^ = AW\m), (G.8) 



with W{t) = U^{t)WUoit). Formal integration leads to 

Jo 

Solving the integral by iteration we obtain 



(G.9) 



1^^(0)1 = 

l + (-zA) / dTW(T) + (-zA)2 / dr C dT'W{T') + . . A \xp{0)). 
^ Jo Jo Jo J 



(G.IO) 



Comparing (G.7 1 and ( G.10| l we obtain the Born expansion: 

U^{t)U;,{t) = 1 + /'dTW(T) 

JO 

+ I dr rdT'W(T)W(T') + ■ ■ ■ . (G.ll) 

JO JO 
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